<!-- HTML header for doxygen 1.8.9.1-->
<!DOCTYPE html PUBLIC "-//W3C//DTD XHTML 1.0 Transitional//EN" "http://www.w3.org/TR/xhtml1/DTD/xhtml1-transitional.dtd">
<html xmlns="http://www.w3.org/1999/xhtml">
<head>
<meta http-equiv="Content-Type" content="text/xhtml;charset=UTF-8"/>
<meta http-equiv="X-UA-Compatible" content="IE=9"/>
<meta name="generator" content="Doxygen 1.8.6"/>
<title>ROSS: core/rand-clcg4.c Source File</title>
<link href="tabs.css" rel="stylesheet" type="text/css"/>
<script type="text/javascript" src="jquery.js"></script>
<script type="text/javascript" src="dynsections.js"></script>
<link href="doxygen.css" rel="stylesheet" type="text/css" />
</head>
<body>
<div id="top"><!-- do not remove this div, it is closed by doxygen! -->
<div id="titlearea">
<table cellspacing="0" cellpadding="0">
 <tbody>
 <tr style="height: 56px;">
  <td style="padding-left: 0.5em;">
   <a href="https://github.com/ROSS-org/ROSS">
   <div id="projectname">ROSS
   </div>
   </a>
  </td>
 </tr>
 </tbody>
</table>
</div>
<!-- end header part -->
<!-- Generated by Doxygen 1.8.6 -->
  <div id="navrow1" class="tabs">
    <ul class="tablist">
      <li><a href="index.html"><span>Main&#160;Page</span></a></li>
      <li><a href="pages.html"><span>Related&#160;Pages</span></a></li>
      <li><a href="annotated.html"><span>Data&#160;Structures</span></a></li>
      <li class="current"><a href="files.html"><span>Files</span></a></li>
    </ul>
  </div>
  <div id="navrow2" class="tabs2">
    <ul class="tablist">
      <li><a href="files.html"><span>File&#160;List</span></a></li>
      <li><a href="globals.html"><span>Globals</span></a></li>
    </ul>
  </div>
<div id="nav-path" class="navpath">
  <ul>
<li class="navelem"><a class="el" href="dir_af18cf245e3394ab3b436135935cfd6e.html">core</a></li>  </ul>
</div>
</div><!-- top -->
<div class="header">
  <div class="headertitle">
<div class="title">rand-clcg4.c</div>  </div>
</div><!--header-->
<div class="contents">
<a href="rand-clcg4_8c.html">Go to the documentation of this file.</a><div class="fragment"><div class="line"><a name="l00001"></a><span class="lineno">    1</span>&#160;<span class="preprocessor">#include &lt;<a class="code" href="ross_8h.html">ross.h</a>&gt;</span></div>
<div class="line"><a name="l00002"></a><span class="lineno">    2</span>&#160;<span class="comment"></span></div>
<div class="line"><a name="l00003"></a><span class="lineno">    3</span>&#160;<span class="comment">/**</span></div>
<div class="line"><a name="l00004"></a><span class="lineno">    4</span>&#160;<span class="comment"> * @file rand-clcg4.c</span></div>
<div class="line"><a name="l00005"></a><span class="lineno">    5</span>&#160;<span class="comment"> * @brief RNG Implementation module</span></div>
<div class="line"><a name="l00006"></a><span class="lineno">    6</span>&#160;<span class="comment"> *</span></div>
<div class="line"><a name="l00007"></a><span class="lineno">    7</span>&#160;<span class="comment"> * Random number generator, provides all of the features GTW requires</span></div>
<div class="line"><a name="l00008"></a><span class="lineno">    8</span>&#160;<span class="comment"> * by default.  Chris hacked this pretty well, he would know the</span></div>
<div class="line"><a name="l00009"></a><span class="lineno">    9</span>&#160;<span class="comment"> * features better.</span></div>
<div class="line"><a name="l00010"></a><span class="lineno">   10</span>&#160;<span class="comment"> *</span></div>
<div class="line"><a name="l00011"></a><span class="lineno">   11</span>&#160;<span class="comment"> */</span></div>
<div class="line"><a name="l00012"></a><span class="lineno">   12</span>&#160;</div>
<div class="line"><a name="l00013"></a><span class="lineno">   13</span>&#160;<span class="comment">/*********************************************************************</span></div>
<div class="line"><a name="l00014"></a><span class="lineno">   14</span>&#160;<span class="comment"> *</span></div>
<div class="line"><a name="l00015"></a><span class="lineno">   15</span>&#160;<span class="comment"> * clcg4 private data/routines</span></div>
<div class="line"><a name="l00016"></a><span class="lineno">   16</span>&#160;<span class="comment"> *</span></div>
<div class="line"><a name="l00017"></a><span class="lineno">   17</span>&#160;<span class="comment"> ********************************************************************/</span></div>
<div class="line"><a name="l00018"></a><span class="lineno">   18</span>&#160;</div>
<div class="line"><a name="l00019"></a><span class="lineno">   19</span>&#160;<span class="comment">// H = 2^15 : use in MultModM.           </span></div>
<div class="line"><a name="l00020"></a><span class="lineno"><a class="line" href="rand-clcg4_8c.html#abec92cc72a096640b821b8cd56a02495">   20</a></span>&#160;<span class="preprocessor">#define H   32768</span></div>
<div class="line"><a name="l00021"></a><span class="lineno">   21</span>&#160;<span class="preprocessor"></span></div>
<div class="line"><a name="l00022"></a><span class="lineno">   22</span>&#160;<span class="comment">// One RNG per PE</span></div>
<div class="line"><a name="l00023"></a><span class="lineno"><a class="line" href="rand-clcg4_8c.html#a5bb18b1128fadea0f468ea648f656ce1">   23</a></span>&#160;<span class="keyword">static</span> <a class="code" href="structtw__rng.html">tw_rng</a>   *<a class="code" href="rand-clcg4_8c.html#a5bb18b1128fadea0f468ea648f656ce1">rng</a> = NULL;</div>
<div class="line"><a name="l00024"></a><span class="lineno">   24</span>&#160;</div>
<div class="line"><a name="l00025"></a><span class="lineno">   25</span>&#160;<span class="comment">// default RNG seed</span></div>
<div class="line"><a name="l00026"></a><span class="lineno"><a class="line" href="rand-clcg4_8c.html#ac6110bfd61d392bab210049c5f6f8fa3">   26</a></span>&#160;int32_t <a class="code" href="rand-clcg4_8c.html#ac6110bfd61d392bab210049c5f6f8fa3">seed</a>[4] = { 11111111, 22222222, 33333333, 44444444 };</div>
<div class="line"><a name="l00027"></a><span class="lineno">   27</span>&#160;<span class="comment"></span></div>
<div class="line"><a name="l00028"></a><span class="lineno">   28</span>&#160;<span class="comment">/**</span></div>
<div class="line"><a name="l00029"></a><span class="lineno">   29</span>&#160;<span class="comment"> * FindB </span></div>
<div class="line"><a name="l00030"></a><span class="lineno">   30</span>&#160;<span class="comment"> * @brief Find B to run CLCG4 backwards</span></div>
<div class="line"><a name="l00031"></a><span class="lineno">   31</span>&#160;<span class="comment"> * </span></div>
<div class="line"><a name="l00032"></a><span class="lineno">   32</span>&#160;<span class="comment"> * B is \f$ a_[i]^{m_[i] - 2} \mathrm{mod m_[i]} \f$                   </span></div>
<div class="line"><a name="l00033"></a><span class="lineno">   33</span>&#160;<span class="comment"> * which is used in running the CLCG4 backwards                        </span></div>
<div class="line"><a name="l00034"></a><span class="lineno">   34</span>&#160;<span class="comment"> * Added by Chris Carothers, 5/15/98                          </span></div>
<div class="line"><a name="l00035"></a><span class="lineno">   35</span>&#160;<span class="comment"> */</span></div>
<div class="line"><a name="l00036"></a><span class="lineno">   36</span>&#160;<span class="keywordtype">long</span> <span class="keywordtype">long</span></div>
<div class="line"><a name="l00037"></a><span class="lineno"><a class="line" href="rand-clcg4_8c.html#ac243c54261ecf17a1b848e9b9a2025a9">   37</a></span>&#160;<a class="code" href="rand-clcg4_8c.html#ac243c54261ecf17a1b848e9b9a2025a9">FindB</a>(<span class="keywordtype">long</span> <span class="keywordtype">long</span> a, <span class="keywordtype">long</span> <span class="keywordtype">long</span> k, <span class="keywordtype">long</span> <span class="keywordtype">long</span> m)</div>
<div class="line"><a name="l00038"></a><span class="lineno">   38</span>&#160;{</div>
<div class="line"><a name="l00039"></a><span class="lineno">   39</span>&#160;  <span class="keywordtype">long</span> <span class="keywordtype">long</span> sqrs[32];</div>
<div class="line"><a name="l00040"></a><span class="lineno">   40</span>&#160;  <span class="keywordtype">long</span> <span class="keywordtype">long</span> power_of_2;</div>
<div class="line"><a name="l00041"></a><span class="lineno">   41</span>&#160;  <span class="keywordtype">long</span> <span class="keywordtype">long</span> b;</div>
<div class="line"><a name="l00042"></a><span class="lineno">   42</span>&#160;</div>
<div class="line"><a name="l00043"></a><span class="lineno">   43</span>&#160;  <span class="keywordtype">int</span> i;</div>
<div class="line"><a name="l00044"></a><span class="lineno">   44</span>&#160;</div>
<div class="line"><a name="l00045"></a><span class="lineno">   45</span>&#160;  sqrs[0] = a;</div>
<div class="line"><a name="l00046"></a><span class="lineno">   46</span>&#160;  <span class="keywordflow">for</span>(i = 1; i &lt; 32; i++)</div>
<div class="line"><a name="l00047"></a><span class="lineno">   47</span>&#160;    {</div>
<div class="line"><a name="l00048"></a><span class="lineno">   48</span>&#160;      sqrs[i] =(sqrs[i - 1] * sqrs[i - 1]) % m;</div>
<div class="line"><a name="l00049"></a><span class="lineno">   49</span>&#160;    }</div>
<div class="line"><a name="l00050"></a><span class="lineno">   50</span>&#160;</div>
<div class="line"><a name="l00051"></a><span class="lineno">   51</span>&#160;  power_of_2 = 1;</div>
<div class="line"><a name="l00052"></a><span class="lineno">   52</span>&#160;  b = 1;</div>
<div class="line"><a name="l00053"></a><span class="lineno">   53</span>&#160;  <span class="keywordflow">for</span>(i = 0; i &lt; 32; i++)</div>
<div class="line"><a name="l00054"></a><span class="lineno">   54</span>&#160;    {</div>
<div class="line"><a name="l00055"></a><span class="lineno">   55</span>&#160;      <span class="keywordflow">if</span>(!(power_of_2 &amp; k))</div>
<div class="line"><a name="l00056"></a><span class="lineno">   56</span>&#160;        {</div>
<div class="line"><a name="l00057"></a><span class="lineno">   57</span>&#160;          sqrs[i] = 1;</div>
<div class="line"><a name="l00058"></a><span class="lineno">   58</span>&#160;</div>
<div class="line"><a name="l00059"></a><span class="lineno">   59</span>&#160;        }</div>
<div class="line"><a name="l00060"></a><span class="lineno">   60</span>&#160;      b =(b * sqrs[i]) % m;</div>
<div class="line"><a name="l00061"></a><span class="lineno">   61</span>&#160;      power_of_2 = power_of_2 * 2;</div>
<div class="line"><a name="l00062"></a><span class="lineno">   62</span>&#160;    }</div>
<div class="line"><a name="l00063"></a><span class="lineno">   63</span>&#160;</div>
<div class="line"><a name="l00064"></a><span class="lineno">   64</span>&#160;  <span class="keywordflow">return</span>(b);</div>
<div class="line"><a name="l00065"></a><span class="lineno">   65</span>&#160;}</div>
<div class="line"><a name="l00066"></a><span class="lineno">   66</span>&#160;<span class="comment"></span></div>
<div class="line"><a name="l00067"></a><span class="lineno">   67</span>&#160;<span class="comment">/**</span></div>
<div class="line"><a name="l00068"></a><span class="lineno">   68</span>&#160;<span class="comment"> * MultiModM </span></div>
<div class="line"><a name="l00069"></a><span class="lineno">   69</span>&#160;<span class="comment"> * @brief Returns(s*t) MOD M                                                        </span></div>
<div class="line"><a name="l00070"></a><span class="lineno">   70</span>&#160;<span class="comment"> * </span></div>
<div class="line"><a name="l00071"></a><span class="lineno">   71</span>&#160;<span class="comment"> * See L&#39;Ecuyer and Cote(1991).                                    </span></div>
<div class="line"><a name="l00072"></a><span class="lineno">   72</span>&#160;<span class="comment"> *</span></div>
<div class="line"><a name="l00073"></a><span class="lineno">   73</span>&#160;<span class="comment"> * Returns(s*t) MOD M.  Assumes that -M &lt; s &lt; M and -M &lt; t &lt; M.    </span></div>
<div class="line"><a name="l00074"></a><span class="lineno">   74</span>&#160;<span class="comment"> */</span></div>
<div class="line"><a name="l00075"></a><span class="lineno">   75</span>&#160;</div>
<div class="line"><a name="l00076"></a><span class="lineno">   76</span>&#160;int32_t</div>
<div class="line"><a name="l00077"></a><span class="lineno"><a class="line" href="rand-clcg4_8c.html#ab89ab01c93d79746cff9a7ed763b24d3">   77</a></span>&#160;<a class="code" href="rand-clcg4_8c.html#ab89ab01c93d79746cff9a7ed763b24d3">MultModM</a>(int32_t s, int32_t t, int32_t M)</div>
<div class="line"><a name="l00078"></a><span class="lineno">   78</span>&#160;{</div>
<div class="line"><a name="l00079"></a><span class="lineno">   79</span>&#160;  int32_t R, S0, S1, q, qh, rh, k;</div>
<div class="line"><a name="l00080"></a><span class="lineno">   80</span>&#160;</div>
<div class="line"><a name="l00081"></a><span class="lineno">   81</span>&#160;  <span class="keywordflow">if</span>(s &lt; 0)</div>
<div class="line"><a name="l00082"></a><span class="lineno">   82</span>&#160;    s += M;</div>
<div class="line"><a name="l00083"></a><span class="lineno">   83</span>&#160;  <span class="keywordflow">if</span>(t &lt; 0)</div>
<div class="line"><a name="l00084"></a><span class="lineno">   84</span>&#160;    t += M;</div>
<div class="line"><a name="l00085"></a><span class="lineno">   85</span>&#160;  <span class="keywordflow">if</span>(s &lt; <a class="code" href="rand-clcg4_8c.html#abec92cc72a096640b821b8cd56a02495">H</a>)</div>
<div class="line"><a name="l00086"></a><span class="lineno">   86</span>&#160;    {</div>
<div class="line"><a name="l00087"></a><span class="lineno">   87</span>&#160;      S0 = s;</div>
<div class="line"><a name="l00088"></a><span class="lineno">   88</span>&#160;      R = 0;</div>
<div class="line"><a name="l00089"></a><span class="lineno">   89</span>&#160;    }</div>
<div class="line"><a name="l00090"></a><span class="lineno">   90</span>&#160;  <span class="keywordflow">else</span></div>
<div class="line"><a name="l00091"></a><span class="lineno">   91</span>&#160;    {</div>
<div class="line"><a name="l00092"></a><span class="lineno">   92</span>&#160;      S1 = s / <a class="code" href="rand-clcg4_8c.html#abec92cc72a096640b821b8cd56a02495">H</a>;</div>
<div class="line"><a name="l00093"></a><span class="lineno">   93</span>&#160;      S0 = s - <a class="code" href="rand-clcg4_8c.html#abec92cc72a096640b821b8cd56a02495">H</a> * S1;</div>
<div class="line"><a name="l00094"></a><span class="lineno">   94</span>&#160;      qh = M / <a class="code" href="rand-clcg4_8c.html#abec92cc72a096640b821b8cd56a02495">H</a>;</div>
<div class="line"><a name="l00095"></a><span class="lineno">   95</span>&#160;      rh = M - <a class="code" href="rand-clcg4_8c.html#abec92cc72a096640b821b8cd56a02495">H</a> * qh;</div>
<div class="line"><a name="l00096"></a><span class="lineno">   96</span>&#160;</div>
<div class="line"><a name="l00097"></a><span class="lineno">   97</span>&#160;      <span class="keywordflow">if</span>(S1 &gt;= <a class="code" href="rand-clcg4_8c.html#abec92cc72a096640b821b8cd56a02495">H</a>)</div>
<div class="line"><a name="l00098"></a><span class="lineno">   98</span>&#160;        {</div>
<div class="line"><a name="l00099"></a><span class="lineno">   99</span>&#160;          S1 -= <a class="code" href="rand-clcg4_8c.html#abec92cc72a096640b821b8cd56a02495">H</a>;</div>
<div class="line"><a name="l00100"></a><span class="lineno">  100</span>&#160;          k = t / qh;</div>
<div class="line"><a name="l00101"></a><span class="lineno">  101</span>&#160;          R = <a class="code" href="rand-clcg4_8c.html#abec92cc72a096640b821b8cd56a02495">H</a> *(t - k * qh) - k * rh;</div>
<div class="line"><a name="l00102"></a><span class="lineno">  102</span>&#160;          <span class="keywordflow">while</span>(R &lt; 0)</div>
<div class="line"><a name="l00103"></a><span class="lineno">  103</span>&#160;            R += M;</div>
<div class="line"><a name="l00104"></a><span class="lineno">  104</span>&#160;        }</div>
<div class="line"><a name="l00105"></a><span class="lineno">  105</span>&#160;      <span class="keywordflow">else</span></div>
<div class="line"><a name="l00106"></a><span class="lineno">  106</span>&#160;        R = 0;</div>
<div class="line"><a name="l00107"></a><span class="lineno">  107</span>&#160;</div>
<div class="line"><a name="l00108"></a><span class="lineno">  108</span>&#160;      <span class="keywordflow">if</span>(S1 != 0)</div>
<div class="line"><a name="l00109"></a><span class="lineno">  109</span>&#160;        {</div>
<div class="line"><a name="l00110"></a><span class="lineno">  110</span>&#160;          q = M / S1;</div>
<div class="line"><a name="l00111"></a><span class="lineno">  111</span>&#160;          k = t / q;</div>
<div class="line"><a name="l00112"></a><span class="lineno">  112</span>&#160;          R -= k *(M - S1 * q);</div>
<div class="line"><a name="l00113"></a><span class="lineno">  113</span>&#160;          <span class="keywordflow">if</span>(R &gt; 0)</div>
<div class="line"><a name="l00114"></a><span class="lineno">  114</span>&#160;            R -= M;</div>
<div class="line"><a name="l00115"></a><span class="lineno">  115</span>&#160;          R += S1 *(t - k * q);</div>
<div class="line"><a name="l00116"></a><span class="lineno">  116</span>&#160;          <span class="keywordflow">while</span>(R &lt; 0)</div>
<div class="line"><a name="l00117"></a><span class="lineno">  117</span>&#160;            R += M;</div>
<div class="line"><a name="l00118"></a><span class="lineno">  118</span>&#160;        }</div>
<div class="line"><a name="l00119"></a><span class="lineno">  119</span>&#160;      k = R / qh;</div>
<div class="line"><a name="l00120"></a><span class="lineno">  120</span>&#160;      R = <a class="code" href="rand-clcg4_8c.html#abec92cc72a096640b821b8cd56a02495">H</a> *(R - k * qh) - k * rh;</div>
<div class="line"><a name="l00121"></a><span class="lineno">  121</span>&#160;      <span class="keywordflow">while</span>(R &lt; 0)</div>
<div class="line"><a name="l00122"></a><span class="lineno">  122</span>&#160;        R += M;</div>
<div class="line"><a name="l00123"></a><span class="lineno">  123</span>&#160;    }</div>
<div class="line"><a name="l00124"></a><span class="lineno">  124</span>&#160;  <span class="keywordflow">if</span>(S0 != 0)</div>
<div class="line"><a name="l00125"></a><span class="lineno">  125</span>&#160;    {</div>
<div class="line"><a name="l00126"></a><span class="lineno">  126</span>&#160;      q = M / S0;</div>
<div class="line"><a name="l00127"></a><span class="lineno">  127</span>&#160;      k = t / q;</div>
<div class="line"><a name="l00128"></a><span class="lineno">  128</span>&#160;      R -= k *(M - S0 * q);</div>
<div class="line"><a name="l00129"></a><span class="lineno">  129</span>&#160;      <span class="keywordflow">if</span>(R &gt; 0)</div>
<div class="line"><a name="l00130"></a><span class="lineno">  130</span>&#160;        R -= M;</div>
<div class="line"><a name="l00131"></a><span class="lineno">  131</span>&#160;      R += S0 *(t - k * q);</div>
<div class="line"><a name="l00132"></a><span class="lineno">  132</span>&#160;      <span class="keywordflow">while</span>(R &lt; 0)</div>
<div class="line"><a name="l00133"></a><span class="lineno">  133</span>&#160;        R += M;</div>
<div class="line"><a name="l00134"></a><span class="lineno">  134</span>&#160;    }</div>
<div class="line"><a name="l00135"></a><span class="lineno">  135</span>&#160;  <span class="keywordflow">return</span> R;</div>
<div class="line"><a name="l00136"></a><span class="lineno">  136</span>&#160;}</div>
<div class="line"><a name="l00137"></a><span class="lineno">  137</span>&#160;</div>
<div class="line"><a name="l00138"></a><span class="lineno">  138</span>&#160;<span class="comment">/********************************************************************</span></div>
<div class="line"><a name="l00139"></a><span class="lineno">  139</span>&#160;<span class="comment"> *</span></div>
<div class="line"><a name="l00140"></a><span class="lineno">  140</span>&#160;<span class="comment"> * clcg4 public interface</span></div>
<div class="line"><a name="l00141"></a><span class="lineno">  141</span>&#160;<span class="comment"> *</span></div>
<div class="line"><a name="l00142"></a><span class="lineno">  142</span>&#160;<span class="comment"> *******************************************************************/</span></div>
<div class="line"><a name="l00143"></a><span class="lineno">  143</span>&#160;</div>
<div class="line"><a name="l00144"></a><span class="lineno"><a class="line" href="rand-clcg4_8c.html#a2888609ddc8a467786e86846bec7ddb9">  144</a></span>&#160;<span class="preprocessor">#define S0_MAX 2147483646</span></div>
<div class="line"><a name="l00145"></a><span class="lineno"><a class="line" href="rand-clcg4_8c.html#a0b3e49b839fa2dd23a2fddf31aaf6ad2">  145</a></span>&#160;<span class="preprocessor"></span><span class="preprocessor">#define S1_MAX 2147483542</span></div>
<div class="line"><a name="l00146"></a><span class="lineno"><a class="line" href="rand-clcg4_8c.html#a5a44a44cb9212f0097f9beec79cc8d41">  146</a></span>&#160;<span class="preprocessor"></span><span class="preprocessor">#define S2_MAX 2147483422</span></div>
<div class="line"><a name="l00147"></a><span class="lineno"><a class="line" href="rand-clcg4_8c.html#ac708ebfb62ad4528bedaeffbe28a3e06">  147</a></span>&#160;<span class="preprocessor"></span><span class="preprocessor">#define S3_MAX 2147483322</span></div>
<div class="line"><a name="l00148"></a><span class="lineno">  148</span>&#160;<span class="preprocessor"></span><span class="comment"></span></div>
<div class="line"><a name="l00149"></a><span class="lineno">  149</span>&#160;<span class="comment">/**</span></div>
<div class="line"><a name="l00150"></a><span class="lineno">  150</span>&#160;<span class="comment"> * These seeds MUST adhere to these requirements.</span></div>
<div class="line"><a name="l00151"></a><span class="lineno">  151</span>&#160;<span class="comment"> * This is explicitly stated in L&#39;Ecuyer and Andres (1997)</span></div>
<div class="line"><a name="l00152"></a><span class="lineno">  152</span>&#160;<span class="comment"> */</span></div>
<div class="line"><a name="l00153"></a><span class="lineno"><a class="line" href="rand-clcg4_8c.html#af6fe13cfedc4e727feee7b28869e036e">  153</a></span>&#160;<span class="keywordtype">void</span> <a class="code" href="rand-clcg4_8c.html#af6fe13cfedc4e727feee7b28869e036e">clamp_seed</a>(uint32_t s[4])</div>
<div class="line"><a name="l00154"></a><span class="lineno">  154</span>&#160;{</div>
<div class="line"><a name="l00155"></a><span class="lineno">  155</span>&#160;    <span class="keywordflow">if</span> (s[0] &lt; 1)</div>
<div class="line"><a name="l00156"></a><span class="lineno">  156</span>&#160;        s[0] = 1;</div>
<div class="line"><a name="l00157"></a><span class="lineno">  157</span>&#160;    <span class="keywordflow">if</span> (s[0] &gt; <a class="code" href="rand-clcg4_8c.html#a2888609ddc8a467786e86846bec7ddb9">S0_MAX</a>)</div>
<div class="line"><a name="l00158"></a><span class="lineno">  158</span>&#160;        s[0] = <a class="code" href="rand-clcg4_8c.html#a2888609ddc8a467786e86846bec7ddb9">S0_MAX</a>;</div>
<div class="line"><a name="l00159"></a><span class="lineno">  159</span>&#160;</div>
<div class="line"><a name="l00160"></a><span class="lineno">  160</span>&#160;    <span class="keywordflow">if</span> (s[1] &lt; 1)</div>
<div class="line"><a name="l00161"></a><span class="lineno">  161</span>&#160;        s[1] = 1;</div>
<div class="line"><a name="l00162"></a><span class="lineno">  162</span>&#160;    <span class="keywordflow">if</span> (s[1] &gt; <a class="code" href="rand-clcg4_8c.html#a0b3e49b839fa2dd23a2fddf31aaf6ad2">S1_MAX</a>)</div>
<div class="line"><a name="l00163"></a><span class="lineno">  163</span>&#160;        s[1] = <a class="code" href="rand-clcg4_8c.html#a0b3e49b839fa2dd23a2fddf31aaf6ad2">S1_MAX</a>;</div>
<div class="line"><a name="l00164"></a><span class="lineno">  164</span>&#160;</div>
<div class="line"><a name="l00165"></a><span class="lineno">  165</span>&#160;    <span class="keywordflow">if</span> (s[2] &lt; 1)</div>
<div class="line"><a name="l00166"></a><span class="lineno">  166</span>&#160;        s[2] = 1;</div>
<div class="line"><a name="l00167"></a><span class="lineno">  167</span>&#160;    <span class="keywordflow">if</span> (s[2] &gt; <a class="code" href="rand-clcg4_8c.html#a5a44a44cb9212f0097f9beec79cc8d41">S2_MAX</a>)</div>
<div class="line"><a name="l00168"></a><span class="lineno">  168</span>&#160;        s[2] = <a class="code" href="rand-clcg4_8c.html#a5a44a44cb9212f0097f9beec79cc8d41">S2_MAX</a>;</div>
<div class="line"><a name="l00169"></a><span class="lineno">  169</span>&#160;</div>
<div class="line"><a name="l00170"></a><span class="lineno">  170</span>&#160;    <span class="keywordflow">if</span> (s[3] &lt; 1)</div>
<div class="line"><a name="l00171"></a><span class="lineno">  171</span>&#160;        s[3] = 1;</div>
<div class="line"><a name="l00172"></a><span class="lineno">  172</span>&#160;    <span class="keywordflow">if</span> (s[3] &gt; <a class="code" href="rand-clcg4_8c.html#ac708ebfb62ad4528bedaeffbe28a3e06">S3_MAX</a>)</div>
<div class="line"><a name="l00173"></a><span class="lineno">  173</span>&#160;        s[3] = <a class="code" href="rand-clcg4_8c.html#ac708ebfb62ad4528bedaeffbe28a3e06">S3_MAX</a>;</div>
<div class="line"><a name="l00174"></a><span class="lineno">  174</span>&#160;}</div>
<div class="line"><a name="l00175"></a><span class="lineno">  175</span>&#160;</div>
<div class="line"><a name="l00176"></a><span class="lineno">  176</span>&#160;<span class="comment">/*</span></div>
<div class="line"><a name="l00177"></a><span class="lineno">  177</span>&#160;<span class="comment"> * rng_set_seed</span></div>
<div class="line"><a name="l00178"></a><span class="lineno">  178</span>&#160;<span class="comment"> */</span></div>
<div class="line"><a name="l00179"></a><span class="lineno">  179</span>&#160;</div>
<div class="line"><a name="l00180"></a><span class="lineno">  180</span>&#160;<span class="keywordtype">void</span></div>
<div class="line"><a name="l00181"></a><span class="lineno"><a class="line" href="rand-clcg4_8c.html#a9c251d8e1876abc122b55fe4156ce14c">  181</a></span>&#160;<a class="code" href="rand-clcg4_8c.html#a9c251d8e1876abc122b55fe4156ce14c">rng_set_seed</a>(<a class="code" href="structtw__rng__stream.html">tw_rng_stream</a> * g, uint32_t s[4])</div>
<div class="line"><a name="l00182"></a><span class="lineno">  182</span>&#160;{</div>
<div class="line"><a name="l00183"></a><span class="lineno">  183</span>&#160;        <span class="keywordtype">int</span>     j;</div>
<div class="line"><a name="l00184"></a><span class="lineno">  184</span>&#160;</div>
<div class="line"><a name="l00185"></a><span class="lineno">  185</span>&#160;    <a class="code" href="rand-clcg4_8c.html#af6fe13cfedc4e727feee7b28869e036e">clamp_seed</a>(s);</div>
<div class="line"><a name="l00186"></a><span class="lineno">  186</span>&#160;</div>
<div class="line"><a name="l00187"></a><span class="lineno">  187</span>&#160;        <span class="keywordflow">for</span>(j = 0; j &lt; 4; j++)</div>
<div class="line"><a name="l00188"></a><span class="lineno">  188</span>&#160;                g-&gt;<a class="code" href="structtw__rng__stream.html#aa90b1f32cf686365db9070d700ca0a38">Ig</a>[j] = s[j];</div>
<div class="line"><a name="l00189"></a><span class="lineno">  189</span>&#160;</div>
<div class="line"><a name="l00190"></a><span class="lineno">  190</span>&#160;        <a class="code" href="rand-clcg4_8c.html#a29a7a9b6a079ced3ba8339a66856fa82">rng_init_generator</a>(g, <a class="code" href="rand-clcg4_8h.html#a3481da76773e3f2621af85e5865484a5a36cbdbd62f02fd7688d6dfc56c5b9670">InitialSeed</a>);</div>
<div class="line"><a name="l00191"></a><span class="lineno">  191</span>&#160;}</div>
<div class="line"><a name="l00192"></a><span class="lineno">  192</span>&#160;</div>
<div class="line"><a name="l00193"></a><span class="lineno">  193</span>&#160;<span class="comment">/*</span></div>
<div class="line"><a name="l00194"></a><span class="lineno">  194</span>&#160;<span class="comment"> * rng_write_state                                                       </span></div>
<div class="line"><a name="l00195"></a><span class="lineno">  195</span>&#160;<span class="comment"> */</span></div>
<div class="line"><a name="l00196"></a><span class="lineno">  196</span>&#160;</div>
<div class="line"><a name="l00197"></a><span class="lineno">  197</span>&#160;<span class="keywordtype">void</span></div>
<div class="line"><a name="l00198"></a><span class="lineno"><a class="line" href="rand-clcg4_8h.html#ad0fa533d1469c99034b45cd302a5134c">  198</a></span>&#160;<a class="code" href="rand-clcg4_8c.html#ad0fa533d1469c99034b45cd302a5134c">rng_write_state</a>(<a class="code" href="structtw__rng__stream.html">tw_rng_stream</a> * g, FILE *f)</div>
<div class="line"><a name="l00199"></a><span class="lineno">  199</span>&#160;{</div>
<div class="line"><a name="l00200"></a><span class="lineno">  200</span>&#160;        <span class="keywordtype">int</span>     j;</div>
<div class="line"><a name="l00201"></a><span class="lineno">  201</span>&#160;</div>
<div class="line"><a name="l00202"></a><span class="lineno">  202</span>&#160;        <span class="keywordflow">for</span>(j = 0; j &lt; 4; j++)</div>
<div class="line"><a name="l00203"></a><span class="lineno">  203</span>&#160;          fprintf( f, <span class="stringliteral">&quot;%u &quot;</span>, g-&gt;<a class="code" href="structtw__rng__stream.html#ae50417100717d3cd46dacb1d7a01e6bf">Cg</a>[j]);</div>
<div class="line"><a name="l00204"></a><span class="lineno">  204</span>&#160;        fprintf( f, <span class="stringliteral">&quot;\n&quot;</span>);</div>
<div class="line"><a name="l00205"></a><span class="lineno">  205</span>&#160;}</div>
<div class="line"><a name="l00206"></a><span class="lineno">  206</span>&#160;</div>
<div class="line"><a name="l00207"></a><span class="lineno">  207</span>&#160;<span class="comment">/*</span></div>
<div class="line"><a name="l00208"></a><span class="lineno">  208</span>&#160;<span class="comment"> * rng_get_state                                                         </span></div>
<div class="line"><a name="l00209"></a><span class="lineno">  209</span>&#160;<span class="comment"> */</span></div>
<div class="line"><a name="l00210"></a><span class="lineno">  210</span>&#160;</div>
<div class="line"><a name="l00211"></a><span class="lineno">  211</span>&#160;<span class="keywordtype">void</span></div>
<div class="line"><a name="l00212"></a><span class="lineno"><a class="line" href="rand-clcg4_8c.html#aaf057a5d3e55f2174f5100a2058f37c5">  212</a></span>&#160;<a class="code" href="rand-clcg4_8c.html#aaf057a5d3e55f2174f5100a2058f37c5">rng_get_state</a>(<a class="code" href="structtw__rng__stream.html">tw_rng_stream</a> * g, uint32_t s[4])</div>
<div class="line"><a name="l00213"></a><span class="lineno">  213</span>&#160;{</div>
<div class="line"><a name="l00214"></a><span class="lineno">  214</span>&#160;        <span class="keywordtype">int</span>     j;</div>
<div class="line"><a name="l00215"></a><span class="lineno">  215</span>&#160;</div>
<div class="line"><a name="l00216"></a><span class="lineno">  216</span>&#160;        <span class="keywordflow">for</span>(j = 0; j &lt; 4; j++)</div>
<div class="line"><a name="l00217"></a><span class="lineno">  217</span>&#160;                s[j] = g-&gt;<a class="code" href="structtw__rng__stream.html#ae50417100717d3cd46dacb1d7a01e6bf">Cg</a>[j];</div>
<div class="line"><a name="l00218"></a><span class="lineno">  218</span>&#160;}</div>
<div class="line"><a name="l00219"></a><span class="lineno">  219</span>&#160;</div>
<div class="line"><a name="l00220"></a><span class="lineno">  220</span>&#160;</div>
<div class="line"><a name="l00221"></a><span class="lineno">  221</span>&#160;<span class="comment">/*</span></div>
<div class="line"><a name="l00222"></a><span class="lineno">  222</span>&#160;<span class="comment"> * rng_get_state                                                         </span></div>
<div class="line"><a name="l00223"></a><span class="lineno">  223</span>&#160;<span class="comment"> */</span></div>
<div class="line"><a name="l00224"></a><span class="lineno">  224</span>&#160;</div>
<div class="line"><a name="l00225"></a><span class="lineno">  225</span>&#160;<span class="keywordtype">void</span></div>
<div class="line"><a name="l00226"></a><span class="lineno"><a class="line" href="rand-clcg4_8c.html#aa97731a64b1cf687109f7a5f8f519c68">  226</a></span>&#160;<a class="code" href="rand-clcg4_8c.html#aa97731a64b1cf687109f7a5f8f519c68">rng_put_state</a>(<a class="code" href="structtw__rng__stream.html">tw_rng_stream</a> * g, uint32_t s[4])</div>
<div class="line"><a name="l00227"></a><span class="lineno">  227</span>&#160;{</div>
<div class="line"><a name="l00228"></a><span class="lineno">  228</span>&#160;        <span class="keywordtype">int</span>     j;</div>
<div class="line"><a name="l00229"></a><span class="lineno">  229</span>&#160;</div>
<div class="line"><a name="l00230"></a><span class="lineno">  230</span>&#160;        <span class="keywordflow">for</span>(j = 0; j &lt; 4; j++)</div>
<div class="line"><a name="l00231"></a><span class="lineno">  231</span>&#160;                g-&gt;<a class="code" href="structtw__rng__stream.html#ae50417100717d3cd46dacb1d7a01e6bf">Cg</a>[j] = s[j];</div>
<div class="line"><a name="l00232"></a><span class="lineno">  232</span>&#160;}</div>
<div class="line"><a name="l00233"></a><span class="lineno">  233</span>&#160;</div>
<div class="line"><a name="l00234"></a><span class="lineno">  234</span>&#160;<span class="comment">/*</span></div>
<div class="line"><a name="l00235"></a><span class="lineno">  235</span>&#160;<span class="comment"> * rng_init_generator                                                    </span></div>
<div class="line"><a name="l00236"></a><span class="lineno">  236</span>&#160;<span class="comment"> */</span></div>
<div class="line"><a name="l00237"></a><span class="lineno">  237</span>&#160;</div>
<div class="line"><a name="l00238"></a><span class="lineno">  238</span>&#160;<span class="keywordtype">void</span></div>
<div class="line"><a name="l00239"></a><span class="lineno"><a class="line" href="rand-clcg4_8h.html#a29a7a9b6a079ced3ba8339a66856fa82">  239</a></span>&#160;<a class="code" href="rand-clcg4_8c.html#a29a7a9b6a079ced3ba8339a66856fa82">rng_init_generator</a>(<a class="code" href="structtw__rng__stream.html">tw_rng_stream</a> * g, <a class="code" href="rand-clcg4_8h.html#a3481da76773e3f2621af85e5865484a5">SeedType</a> Where)</div>
<div class="line"><a name="l00240"></a><span class="lineno">  240</span>&#160;{</div>
<div class="line"><a name="l00241"></a><span class="lineno">  241</span>&#160;        <span class="keywordtype">int</span>     j;</div>
<div class="line"><a name="l00242"></a><span class="lineno">  242</span>&#160;</div>
<div class="line"><a name="l00243"></a><span class="lineno">  243</span>&#160;        <span class="keywordflow">for</span>(j = 0; j &lt; 4; j++)</div>
<div class="line"><a name="l00244"></a><span class="lineno">  244</span>&#160;        {</div>
<div class="line"><a name="l00245"></a><span class="lineno">  245</span>&#160;                <span class="keywordflow">switch</span>(Where)</div>
<div class="line"><a name="l00246"></a><span class="lineno">  246</span>&#160;                {</div>
<div class="line"><a name="l00247"></a><span class="lineno">  247</span>&#160;                        <span class="keywordflow">case</span> <a class="code" href="rand-clcg4_8h.html#a3481da76773e3f2621af85e5865484a5a36cbdbd62f02fd7688d6dfc56c5b9670">InitialSeed</a>:</div>
<div class="line"><a name="l00248"></a><span class="lineno">  248</span>&#160;                                g-&gt;<a class="code" href="structtw__rng__stream.html#a61563015a585f69852a131cb62a4a5b1">Lg</a>[j] = g-&gt;<a class="code" href="structtw__rng__stream.html#aa90b1f32cf686365db9070d700ca0a38">Ig</a>[j];</div>
<div class="line"><a name="l00249"></a><span class="lineno">  249</span>&#160;                                <span class="keywordflow">break</span>;</div>
<div class="line"><a name="l00250"></a><span class="lineno">  250</span>&#160;                        <span class="keywordflow">case</span> <a class="code" href="rand-clcg4_8h.html#a3481da76773e3f2621af85e5865484a5a6f0f1982e69fd30c91d7d3b7097f8912">NewSeed</a>:</div>
<div class="line"><a name="l00251"></a><span class="lineno">  251</span>&#160;                                g-&gt;<a class="code" href="structtw__rng__stream.html#a61563015a585f69852a131cb62a4a5b1">Lg</a>[j] = <a class="code" href="rand-clcg4_8c.html#ab89ab01c93d79746cff9a7ed763b24d3">MultModM</a>(rng-&gt;<a class="code" href="structtw__rng.html#ad64150f8992ffc7cff059682f5418378">aw</a>[j], g-&gt;<a class="code" href="structtw__rng__stream.html#a61563015a585f69852a131cb62a4a5b1">Lg</a>[j], rng-&gt;<a class="code" href="structtw__rng.html#ab6cd71dbaaf9c25885b5647c4e1f9d02">m</a>[j]);</div>
<div class="line"><a name="l00252"></a><span class="lineno">  252</span>&#160;                                <span class="keywordflow">break</span>;</div>
<div class="line"><a name="l00253"></a><span class="lineno">  253</span>&#160;                        <span class="keywordflow">case</span> <a class="code" href="rand-clcg4_8h.html#a3481da76773e3f2621af85e5865484a5a87c13b1e6a443e099ca4aa828db9e7a3">LastSeed</a>:</div>
<div class="line"><a name="l00254"></a><span class="lineno">  254</span>&#160;                                <span class="keywordflow">break</span>;</div>
<div class="line"><a name="l00255"></a><span class="lineno">  255</span>&#160;                }</div>
<div class="line"><a name="l00256"></a><span class="lineno">  256</span>&#160;</div>
<div class="line"><a name="l00257"></a><span class="lineno">  257</span>&#160;                g-&gt;<a class="code" href="structtw__rng__stream.html#ae50417100717d3cd46dacb1d7a01e6bf">Cg</a>[j] = g-&gt;<a class="code" href="structtw__rng__stream.html#a61563015a585f69852a131cb62a4a5b1">Lg</a>[j];</div>
<div class="line"><a name="l00258"></a><span class="lineno">  258</span>&#160;        }</div>
<div class="line"><a name="l00259"></a><span class="lineno">  259</span>&#160;}</div>
<div class="line"><a name="l00260"></a><span class="lineno">  260</span>&#160;</div>
<div class="line"><a name="l00261"></a><span class="lineno">  261</span>&#160;<span class="comment">/*</span></div>
<div class="line"><a name="l00262"></a><span class="lineno">  262</span>&#160;<span class="comment"> * rng_set_initial_seed                                                   </span></div>
<div class="line"><a name="l00263"></a><span class="lineno">  263</span>&#160;<span class="comment"> */</span></div>
<div class="line"><a name="l00264"></a><span class="lineno">  264</span>&#160;<span class="keywordtype">void</span></div>
<div class="line"><a name="l00265"></a><span class="lineno"><a class="line" href="ross-random_8h.html#a9179f0148c620729309111ad61a3fa5b">  265</a></span>&#160;<a class="code" href="rand-clcg4_8c.html#a9179f0148c620729309111ad61a3fa5b">tw_rand_initial_seed</a>(<a class="code" href="structtw__rng__stream.html">tw_rng_stream</a> * g, <a class="code" href="ross_8h.html#a911515dff2fb79886a16ed44df5f9d20">tw_lpid</a> <span class="keywordtype">id</span>)</div>
<div class="line"><a name="l00266"></a><span class="lineno">  266</span>&#160;{</div>
<div class="line"><a name="l00267"></a><span class="lineno">  267</span>&#160;        <a class="code" href="ross_8h.html#a911515dff2fb79886a16ed44df5f9d20">tw_lpid</a> mask_bit = 1;</div>
<div class="line"><a name="l00268"></a><span class="lineno">  268</span>&#160;</div>
<div class="line"><a name="l00269"></a><span class="lineno">  269</span>&#160;        uint32_t Ig_t[4];</div>
<div class="line"><a name="l00270"></a><span class="lineno">  270</span>&#160;        uint32_t avw_t[4];</div>
<div class="line"><a name="l00271"></a><span class="lineno">  271</span>&#160;</div>
<div class="line"><a name="l00272"></a><span class="lineno">  272</span>&#160;        <span class="keywordtype">int</span> i;</div>
<div class="line"><a name="l00273"></a><span class="lineno">  273</span>&#160;        <span class="keywordtype">int</span> j;</div>
<div class="line"><a name="l00274"></a><span class="lineno">  274</span>&#160;        <span class="keywordtype">int</span> positions = ((<span class="keyword">sizeof</span>(<a class="code" href="ross_8h.html#a911515dff2fb79886a16ed44df5f9d20">tw_lpid</a>)) * 8) - 1;</div>
<div class="line"><a name="l00275"></a><span class="lineno">  275</span>&#160;</div>
<div class="line"><a name="l00276"></a><span class="lineno">  276</span>&#160;        <span class="comment">//seed for zero</span></div>
<div class="line"><a name="l00277"></a><span class="lineno">  277</span>&#160;        <span class="keywordflow">for</span>(j = 0; j &lt; 4; j++)</div>
<div class="line"><a name="l00278"></a><span class="lineno">  278</span>&#160;                Ig_t[j] = rng-&gt;<a class="code" href="structtw__rng.html#ada588d3e377c9f0d10fb22406f5b47c6">seed</a>[j];</div>
<div class="line"><a name="l00279"></a><span class="lineno">  279</span>&#160;</div>
<div class="line"><a name="l00280"></a><span class="lineno">  280</span>&#160;        mask_bit &lt;&lt;= positions;</div>
<div class="line"><a name="l00281"></a><span class="lineno">  281</span>&#160;</div>
<div class="line"><a name="l00282"></a><span class="lineno">  282</span>&#160;        g-&gt;<a class="code" href="structtw__rng__stream.html#ae373cfc091d1ac7361776979d914ad59">count</a> = 0;</div>
<div class="line"><a name="l00283"></a><span class="lineno">  283</span>&#160;</div>
<div class="line"><a name="l00284"></a><span class="lineno">  284</span>&#160;        <span class="keywordflow">do</span></div>
<div class="line"><a name="l00285"></a><span class="lineno">  285</span>&#160;        {</div>
<div class="line"><a name="l00286"></a><span class="lineno">  286</span>&#160;                <span class="keywordflow">if</span>(<span class="keywordtype">id</span> &amp; mask_bit)</div>
<div class="line"><a name="l00287"></a><span class="lineno">  287</span>&#160;                {</div>
<div class="line"><a name="l00288"></a><span class="lineno">  288</span>&#160;                        <span class="keywordflow">for</span>(j = 0; j &lt; 4; j++)</div>
<div class="line"><a name="l00289"></a><span class="lineno">  289</span>&#160;                        {</div>
<div class="line"><a name="l00290"></a><span class="lineno">  290</span>&#160;                                avw_t[j] = rng-&gt;<a class="code" href="structtw__rng.html#a4f5b2d7e79fcebe5c9731b1cf3ee98db">avw</a>[j];</div>
<div class="line"><a name="l00291"></a><span class="lineno">  291</span>&#160;</div>
<div class="line"><a name="l00292"></a><span class="lineno">  292</span>&#160;                                <span class="comment">// exponentiate modulus</span></div>
<div class="line"><a name="l00293"></a><span class="lineno">  293</span>&#160;                                <span class="keywordflow">for</span>(i = 0; i &lt; positions; i++)</div>
<div class="line"><a name="l00294"></a><span class="lineno">  294</span>&#160;                                        avw_t[j] = <a class="code" href="rand-clcg4_8c.html#ab89ab01c93d79746cff9a7ed763b24d3">MultModM</a>(avw_t[j], avw_t[j], rng-&gt;<a class="code" href="structtw__rng.html#ab6cd71dbaaf9c25885b5647c4e1f9d02">m</a>[j]);</div>
<div class="line"><a name="l00295"></a><span class="lineno">  295</span>&#160;</div>
<div class="line"><a name="l00296"></a><span class="lineno">  296</span>&#160;                                Ig_t[j] = <a class="code" href="rand-clcg4_8c.html#ab89ab01c93d79746cff9a7ed763b24d3">MultModM</a>(avw_t[j], Ig_t[j], rng-&gt;<a class="code" href="structtw__rng.html#ab6cd71dbaaf9c25885b5647c4e1f9d02">m</a>[j]);</div>
<div class="line"><a name="l00297"></a><span class="lineno">  297</span>&#160;                        }</div>
<div class="line"><a name="l00298"></a><span class="lineno">  298</span>&#160;                }</div>
<div class="line"><a name="l00299"></a><span class="lineno">  299</span>&#160;</div>
<div class="line"><a name="l00300"></a><span class="lineno">  300</span>&#160;                mask_bit &gt;&gt;= 1;</div>
<div class="line"><a name="l00301"></a><span class="lineno">  301</span>&#160;                positions--;</div>
<div class="line"><a name="l00302"></a><span class="lineno">  302</span>&#160;        } <span class="keywordflow">while</span>(positions &gt; 0);</div>
<div class="line"><a name="l00303"></a><span class="lineno">  303</span>&#160;</div>
<div class="line"><a name="l00304"></a><span class="lineno">  304</span>&#160;        <span class="keywordflow">if</span>(<span class="keywordtype">id</span> % 2)</div>
<div class="line"><a name="l00305"></a><span class="lineno">  305</span>&#160;        {</div>
<div class="line"><a name="l00306"></a><span class="lineno">  306</span>&#160;                <span class="keywordflow">for</span>(j = 0; j &lt; 4; j++)</div>
<div class="line"><a name="l00307"></a><span class="lineno">  307</span>&#160;                        Ig_t[j] = <a class="code" href="rand-clcg4_8c.html#ab89ab01c93d79746cff9a7ed763b24d3">MultModM</a>(rng-&gt;<a class="code" href="structtw__rng.html#a4f5b2d7e79fcebe5c9731b1cf3ee98db">avw</a>[j], Ig_t[j], rng-&gt;<a class="code" href="structtw__rng.html#ab6cd71dbaaf9c25885b5647c4e1f9d02">m</a>[j]);</div>
<div class="line"><a name="l00308"></a><span class="lineno">  308</span>&#160;        }</div>
<div class="line"><a name="l00309"></a><span class="lineno">  309</span>&#160;</div>
<div class="line"><a name="l00310"></a><span class="lineno">  310</span>&#160;    <a class="code" href="rand-clcg4_8c.html#af6fe13cfedc4e727feee7b28869e036e">clamp_seed</a>(Ig_t);</div>
<div class="line"><a name="l00311"></a><span class="lineno">  311</span>&#160;</div>
<div class="line"><a name="l00312"></a><span class="lineno">  312</span>&#160;        <span class="keywordflow">for</span>(j = 0; j &lt; 4; j++)</div>
<div class="line"><a name="l00313"></a><span class="lineno">  313</span>&#160;                g-&gt;Ig[j] = Ig_t[j];</div>
<div class="line"><a name="l00314"></a><span class="lineno">  314</span>&#160;</div>
<div class="line"><a name="l00315"></a><span class="lineno">  315</span>&#160;        <a class="code" href="rand-clcg4_8c.html#a29a7a9b6a079ced3ba8339a66856fa82">rng_init_generator</a>(g, <a class="code" href="rand-clcg4_8h.html#a3481da76773e3f2621af85e5865484a5a36cbdbd62f02fd7688d6dfc56c5b9670">InitialSeed</a>);</div>
<div class="line"><a name="l00316"></a><span class="lineno">  316</span>&#160;        <span class="comment">//rng_write_state(g);</span></div>
<div class="line"><a name="l00317"></a><span class="lineno">  317</span>&#160;}</div>
<div class="line"><a name="l00318"></a><span class="lineno">  318</span>&#160;</div>
<div class="line"><a name="l00319"></a><span class="lineno">  319</span>&#160;<span class="keywordtype">void</span></div>
<div class="line"><a name="l00320"></a><span class="lineno"><a class="line" href="ross-extern_8h.html#a96c2fb54945c414a31ce0512c39622b0">  320</a></span>&#160;<a class="code" href="rand-clcg4_8c.html#a96c2fb54945c414a31ce0512c39622b0">tw_rand_init_streams</a>(<a class="code" href="structtw__lp.html">tw_lp</a> * lp, <span class="keywordtype">unsigned</span> <span class="keywordtype">int</span> nstreams)</div>
<div class="line"><a name="l00321"></a><span class="lineno">  321</span>&#160;{</div>
<div class="line"><a name="l00322"></a><span class="lineno">  322</span>&#160;        <span class="keywordtype">unsigned</span> <span class="keywordtype">int</span> i;</div>
<div class="line"><a name="l00323"></a><span class="lineno">  323</span>&#160;</div>
<div class="line"><a name="l00324"></a><span class="lineno">  324</span>&#160;        <span class="keywordflow">if</span>(nstreams &gt; <a class="code" href="ross-extern_8h.html#a99999f4e30648f5007c61723c95a2d0e">g_tw_nRNG_per_lp</a>)</div>
<div class="line"><a name="l00325"></a><span class="lineno">  325</span>&#160;            <a class="code" href="ross-extern_8h.html#a49ed2388aaae26e43280b7909c834aaa">tw_error</a>(<a class="code" href="ross-extern_8h.html#afffdf140f2c492f61546ad783833f127">TW_LOC</a>, <span class="stringliteral">&quot;LP %lu asked for more RNG streams (%d) than the global maximum (g_tw_nRNG_per_lp:%d)\n&quot;</span>, lp-&gt;<a class="code" href="structtw__lp.html#a3e75d9bf7796c16870ad686564ee4c20">gid</a>, nstreams, <a class="code" href="ross-extern_8h.html#a99999f4e30648f5007c61723c95a2d0e">g_tw_nRNG_per_lp</a>);</div>
<div class="line"><a name="l00326"></a><span class="lineno">  326</span>&#160;</div>
<div class="line"><a name="l00327"></a><span class="lineno">  327</span>&#160;        lp-&gt;<a class="code" href="structtw__lp.html#af7e48c8a55f5fff35746ac044df1902f">rng</a> = (<a class="code" href="structtw__rng__stream.html">tw_rng_stream</a> *) <a class="code" href="ross-extern_8h.html#a445b1ed30cd2e8ef1f837e782a339a70">tw_calloc</a>(<a class="code" href="ross-extern_8h.html#afffdf140f2c492f61546ad783833f127">TW_LOC</a>, <span class="stringliteral">&quot;LP RNG Streams&quot;</span>, <span class="keyword">sizeof</span>(*lp-&gt;<a class="code" href="structtw__lp.html#af7e48c8a55f5fff35746ac044df1902f">rng</a>), nstreams);</div>
<div class="line"><a name="l00328"></a><span class="lineno">  328</span>&#160;</div>
<div class="line"><a name="l00329"></a><span class="lineno">  329</span>&#160;        <span class="keywordflow">for</span>(i = 0; i &lt; nstreams; i++) {</div>
<div class="line"><a name="l00330"></a><span class="lineno">  330</span>&#160;            <a class="code" href="rand-clcg4_8c.html#a9179f0148c620729309111ad61a3fa5b">tw_rand_initial_seed</a>(&amp;lp-&gt;<a class="code" href="structtw__lp.html#af7e48c8a55f5fff35746ac044df1902f">rng</a>[i], (lp-&gt;<a class="code" href="structtw__lp.html#a3e75d9bf7796c16870ad686564ee4c20">gid</a> * <a class="code" href="ross-extern_8h.html#a99999f4e30648f5007c61723c95a2d0e">g_tw_nRNG_per_lp</a>) + i);</div>
<div class="line"><a name="l00331"></a><span class="lineno">  331</span>&#160;        }</div>
<div class="line"><a name="l00332"></a><span class="lineno">  332</span>&#160;}</div>
<div class="line"><a name="l00333"></a><span class="lineno">  333</span>&#160;</div>
<div class="line"><a name="l00334"></a><span class="lineno">  334</span>&#160;<span class="comment">/*</span></div>
<div class="line"><a name="l00335"></a><span class="lineno">  335</span>&#160;<span class="comment"> * rng_init                                                             </span></div>
<div class="line"><a name="l00336"></a><span class="lineno">  336</span>&#160;<span class="comment"> */</span></div>
<div class="line"><a name="l00337"></a><span class="lineno">  337</span>&#160;<a class="code" href="structtw__rng.html">tw_rng</a>  *</div>
<div class="line"><a name="l00338"></a><span class="lineno"><a class="line" href="rand-clcg4_8h.html#a3c5ea1f9125b9534b62a80dbfa20dafb">  338</a></span>&#160;<a class="code" href="rand-clcg4_8c.html#a3c5ea1f9125b9534b62a80dbfa20dafb">rng_init</a>(<span class="keywordtype">int</span> v, <span class="keywordtype">int</span> w)</div>
<div class="line"><a name="l00339"></a><span class="lineno">  339</span>&#160;{</div>
<div class="line"><a name="l00340"></a><span class="lineno">  340</span>&#160;        <span class="keywordtype">int</span>      i;</div>
<div class="line"><a name="l00341"></a><span class="lineno">  341</span>&#160;        <span class="keywordtype">int</span>      j;</div>
<div class="line"><a name="l00342"></a><span class="lineno">  342</span>&#160;</div>
<div class="line"><a name="l00343"></a><span class="lineno">  343</span>&#160;        rng = (<a class="code" href="structtw__rng.html">tw_rng</a> *) <a class="code" href="ross-extern_8h.html#a445b1ed30cd2e8ef1f837e782a339a70">tw_calloc</a>(<a class="code" href="ross-extern_8h.html#afffdf140f2c492f61546ad783833f127">TW_LOC</a>, <span class="stringliteral">&quot;RNG&quot;</span>, <span class="keyword">sizeof</span>(*rng), 1);</div>
<div class="line"><a name="l00344"></a><span class="lineno">  344</span>&#160;</div>
<div class="line"><a name="l00345"></a><span class="lineno">  345</span>&#160;        rng-&gt;<a class="code" href="structtw__rng.html#ab6cd71dbaaf9c25885b5647c4e1f9d02">m</a>[0] = 2147483647;</div>
<div class="line"><a name="l00346"></a><span class="lineno">  346</span>&#160;        rng-&gt;<a class="code" href="structtw__rng.html#ab6cd71dbaaf9c25885b5647c4e1f9d02">m</a>[1] = 2147483543;</div>
<div class="line"><a name="l00347"></a><span class="lineno">  347</span>&#160;        rng-&gt;<a class="code" href="structtw__rng.html#ab6cd71dbaaf9c25885b5647c4e1f9d02">m</a>[2] = 2147483423;</div>
<div class="line"><a name="l00348"></a><span class="lineno">  348</span>&#160;        rng-&gt;<a class="code" href="structtw__rng.html#ab6cd71dbaaf9c25885b5647c4e1f9d02">m</a>[3] = 2147483323;</div>
<div class="line"><a name="l00349"></a><span class="lineno">  349</span>&#160;</div>
<div class="line"><a name="l00350"></a><span class="lineno">  350</span>&#160;        rng-&gt;<a class="code" href="structtw__rng.html#aff3f20edc0b90fd3191487c21b2de697">a</a>[0] = 45991;</div>
<div class="line"><a name="l00351"></a><span class="lineno">  351</span>&#160;        rng-&gt;<a class="code" href="structtw__rng.html#aff3f20edc0b90fd3191487c21b2de697">a</a>[1] = 207707;</div>
<div class="line"><a name="l00352"></a><span class="lineno">  352</span>&#160;        rng-&gt;<a class="code" href="structtw__rng.html#aff3f20edc0b90fd3191487c21b2de697">a</a>[2] = 138556;</div>
<div class="line"><a name="l00353"></a><span class="lineno">  353</span>&#160;        rng-&gt;<a class="code" href="structtw__rng.html#aff3f20edc0b90fd3191487c21b2de697">a</a>[3] = 49689;</div>
<div class="line"><a name="l00354"></a><span class="lineno">  354</span>&#160;        </div>
<div class="line"><a name="l00355"></a><span class="lineno">  355</span>&#160;        <span class="keywordflow">if</span>(<a class="code" href="ross-extern_8h.html#a11d3d6ea7b82645df5fc7011a407c355">g_tw_rng_seed</a>)</div>
<div class="line"><a name="l00356"></a><span class="lineno">  356</span>&#160;        {</div>
<div class="line"><a name="l00357"></a><span class="lineno">  357</span>&#160;        <a class="code" href="rand-clcg4_8c.html#af6fe13cfedc4e727feee7b28869e036e">clamp_seed</a>((uint32_t*)<a class="code" href="ross-extern_8h.html#a11d3d6ea7b82645df5fc7011a407c355">g_tw_rng_seed</a>);</div>
<div class="line"><a name="l00358"></a><span class="lineno">  358</span>&#160;                <span class="keywordflow">for</span>(j = 0; j &lt; 4; j++)</div>
<div class="line"><a name="l00359"></a><span class="lineno">  359</span>&#160;                        rng-&gt;<a class="code" href="structtw__rng.html#ada588d3e377c9f0d10fb22406f5b47c6">seed</a>[j] = <a class="code" href="ross-extern_8h.html#a11d3d6ea7b82645df5fc7011a407c355">g_tw_rng_seed</a>[j];</div>
<div class="line"><a name="l00360"></a><span class="lineno">  360</span>&#160;        } <span class="keywordflow">else</span></div>
<div class="line"><a name="l00361"></a><span class="lineno">  361</span>&#160;        {</div>
<div class="line"><a name="l00362"></a><span class="lineno">  362</span>&#160;                rng-&gt;<a class="code" href="structtw__rng.html#ada588d3e377c9f0d10fb22406f5b47c6">seed</a>[0] = 11111111;</div>
<div class="line"><a name="l00363"></a><span class="lineno">  363</span>&#160;                rng-&gt;<a class="code" href="structtw__rng.html#ada588d3e377c9f0d10fb22406f5b47c6">seed</a>[1] = 22222222;</div>
<div class="line"><a name="l00364"></a><span class="lineno">  364</span>&#160;                rng-&gt;<a class="code" href="structtw__rng.html#ada588d3e377c9f0d10fb22406f5b47c6">seed</a>[2] = 33333333;</div>
<div class="line"><a name="l00365"></a><span class="lineno">  365</span>&#160;                rng-&gt;<a class="code" href="structtw__rng.html#ada588d3e377c9f0d10fb22406f5b47c6">seed</a>[3] = 44444444;</div>
<div class="line"><a name="l00366"></a><span class="lineno">  366</span>&#160;        }</div>
<div class="line"><a name="l00367"></a><span class="lineno">  367</span>&#160;</div>
<div class="line"><a name="l00368"></a><span class="lineno">  368</span>&#160;        <span class="keywordflow">for</span>(j = 0; j &lt; 4; j++)</div>
<div class="line"><a name="l00369"></a><span class="lineno">  369</span>&#160;                rng-&gt;<a class="code" href="structtw__rng.html#ad64150f8992ffc7cff059682f5418378">aw</a>[j] = rng-&gt;<a class="code" href="structtw__rng.html#aff3f20edc0b90fd3191487c21b2de697">a</a>[j];</div>
<div class="line"><a name="l00370"></a><span class="lineno">  370</span>&#160;</div>
<div class="line"><a name="l00371"></a><span class="lineno">  371</span>&#160;        <span class="keywordflow">for</span>(j = 0; j &lt; 4; j++)</div>
<div class="line"><a name="l00372"></a><span class="lineno">  372</span>&#160;        {</div>
<div class="line"><a name="l00373"></a><span class="lineno">  373</span>&#160;                <span class="keywordflow">for</span>(i = 1; i &lt;= w; i++)</div>
<div class="line"><a name="l00374"></a><span class="lineno">  374</span>&#160;                        rng-&gt;<a class="code" href="structtw__rng.html#ad64150f8992ffc7cff059682f5418378">aw</a>[j] = <a class="code" href="rand-clcg4_8c.html#ab89ab01c93d79746cff9a7ed763b24d3">MultModM</a>(rng-&gt;<a class="code" href="structtw__rng.html#ad64150f8992ffc7cff059682f5418378">aw</a>[j], rng-&gt;<a class="code" href="structtw__rng.html#ad64150f8992ffc7cff059682f5418378">aw</a>[j], rng-&gt;<a class="code" href="structtw__rng.html#ab6cd71dbaaf9c25885b5647c4e1f9d02">m</a>[j]);</div>
<div class="line"><a name="l00375"></a><span class="lineno">  375</span>&#160;</div>
<div class="line"><a name="l00376"></a><span class="lineno">  376</span>&#160;                rng-&gt;<a class="code" href="structtw__rng.html#a4f5b2d7e79fcebe5c9731b1cf3ee98db">avw</a>[j] = rng-&gt;<a class="code" href="structtw__rng.html#ad64150f8992ffc7cff059682f5418378">aw</a>[j];</div>
<div class="line"><a name="l00377"></a><span class="lineno">  377</span>&#160;</div>
<div class="line"><a name="l00378"></a><span class="lineno">  378</span>&#160;                <span class="keywordflow">for</span>(i = 1; i &lt;= v; i++)</div>
<div class="line"><a name="l00379"></a><span class="lineno">  379</span>&#160;                        rng-&gt;<a class="code" href="structtw__rng.html#a4f5b2d7e79fcebe5c9731b1cf3ee98db">avw</a>[j] = <a class="code" href="rand-clcg4_8c.html#ab89ab01c93d79746cff9a7ed763b24d3">MultModM</a>(rng-&gt;<a class="code" href="structtw__rng.html#a4f5b2d7e79fcebe5c9731b1cf3ee98db">avw</a>[j], rng-&gt;<a class="code" href="structtw__rng.html#a4f5b2d7e79fcebe5c9731b1cf3ee98db">avw</a>[j], rng-&gt;<a class="code" href="structtw__rng.html#ab6cd71dbaaf9c25885b5647c4e1f9d02">m</a>[j]);</div>
<div class="line"><a name="l00380"></a><span class="lineno">  380</span>&#160;        }</div>
<div class="line"><a name="l00381"></a><span class="lineno">  381</span>&#160;</div>
<div class="line"><a name="l00382"></a><span class="lineno">  382</span>&#160;        <span class="keywordflow">for</span>(j = 0; j &lt; 4; j++)</div>
<div class="line"><a name="l00383"></a><span class="lineno">  383</span>&#160;                rng-&gt;<a class="code" href="structtw__rng.html#a3966a408c9f22a22d2e7ebc387156086">b</a>[j] = <a class="code" href="rand-clcg4_8c.html#ac243c54261ecf17a1b848e9b9a2025a9">FindB</a>(rng-&gt;<a class="code" href="structtw__rng.html#aff3f20edc0b90fd3191487c21b2de697">a</a>[j],(rng-&gt;<a class="code" href="structtw__rng.html#ab6cd71dbaaf9c25885b5647c4e1f9d02">m</a>[j] - 2), rng-&gt;<a class="code" href="structtw__rng.html#ab6cd71dbaaf9c25885b5647c4e1f9d02">m</a>[j]);</div>
<div class="line"><a name="l00384"></a><span class="lineno">  384</span>&#160;</div>
<div class="line"><a name="l00385"></a><span class="lineno">  385</span>&#160;        <span class="keywordflow">return</span> <a class="code" href="rand-clcg4_8c.html#a5bb18b1128fadea0f468ea648f656ce1">rng</a>;</div>
<div class="line"><a name="l00386"></a><span class="lineno">  386</span>&#160;}</div>
<div class="line"><a name="l00387"></a><span class="lineno">  387</span>&#160;</div>
<div class="line"><a name="l00388"></a><span class="lineno">  388</span>&#160;<span class="comment">/*</span></div>
<div class="line"><a name="l00389"></a><span class="lineno">  389</span>&#160;<span class="comment"> * rng_gen_val                                                              </span></div>
<div class="line"><a name="l00390"></a><span class="lineno">  390</span>&#160;<span class="comment"> */</span></div>
<div class="line"><a name="l00391"></a><span class="lineno">  391</span>&#160;<span class="keywordtype">double</span></div>
<div class="line"><a name="l00392"></a><span class="lineno"><a class="line" href="rand-clcg4_8h.html#adc7b882f7220d319c89b7500c5d2395f">  392</a></span>&#160;<a class="code" href="rand-clcg4_8c.html#adc7b882f7220d319c89b7500c5d2395f">rng_gen_val</a>(<a class="code" href="structtw__rng__stream.html">tw_rng_stream</a> * g)</div>
<div class="line"><a name="l00393"></a><span class="lineno">  393</span>&#160;{</div>
<div class="line"><a name="l00394"></a><span class="lineno">  394</span>&#160;  int32_t k, s;</div>
<div class="line"><a name="l00395"></a><span class="lineno">  395</span>&#160;  <span class="keywordtype">double</span> u;</div>
<div class="line"><a name="l00396"></a><span class="lineno">  396</span>&#160;</div>
<div class="line"><a name="l00397"></a><span class="lineno">  397</span>&#160;  g-&gt;<a class="code" href="structtw__rng__stream.html#ae373cfc091d1ac7361776979d914ad59">count</a>++;</div>
<div class="line"><a name="l00398"></a><span class="lineno">  398</span>&#160;</div>
<div class="line"><a name="l00399"></a><span class="lineno">  399</span>&#160;  u = 0.0;</div>
<div class="line"><a name="l00400"></a><span class="lineno">  400</span>&#160;</div>
<div class="line"><a name="l00401"></a><span class="lineno">  401</span>&#160;  s = g-&gt;<a class="code" href="structtw__rng__stream.html#ae50417100717d3cd46dacb1d7a01e6bf">Cg</a>[0];</div>
<div class="line"><a name="l00402"></a><span class="lineno">  402</span>&#160;  k = s / 46693;</div>
<div class="line"><a name="l00403"></a><span class="lineno">  403</span>&#160;  s = 45991 *(s - k * 46693) - k * 25884;</div>
<div class="line"><a name="l00404"></a><span class="lineno">  404</span>&#160;  <span class="keywordflow">if</span>(s &lt; 0)</div>
<div class="line"><a name="l00405"></a><span class="lineno">  405</span>&#160;    s = s + 2147483647;</div>
<div class="line"><a name="l00406"></a><span class="lineno">  406</span>&#160;  g-&gt;<a class="code" href="structtw__rng__stream.html#ae50417100717d3cd46dacb1d7a01e6bf">Cg</a>[0] = s;</div>
<div class="line"><a name="l00407"></a><span class="lineno">  407</span>&#160;  u = u + 4.65661287524579692e-10 * s;</div>
<div class="line"><a name="l00408"></a><span class="lineno">  408</span>&#160;</div>
<div class="line"><a name="l00409"></a><span class="lineno">  409</span>&#160;  s = g-&gt;<a class="code" href="structtw__rng__stream.html#ae50417100717d3cd46dacb1d7a01e6bf">Cg</a>[1];</div>
<div class="line"><a name="l00410"></a><span class="lineno">  410</span>&#160;  k = s / 10339;</div>
<div class="line"><a name="l00411"></a><span class="lineno">  411</span>&#160;  s = 207707 *(s - k * 10339) - k * 870;</div>
<div class="line"><a name="l00412"></a><span class="lineno">  412</span>&#160;  <span class="keywordflow">if</span>(s &lt; 0)</div>
<div class="line"><a name="l00413"></a><span class="lineno">  413</span>&#160;    s = s + 2147483543;</div>
<div class="line"><a name="l00414"></a><span class="lineno">  414</span>&#160;  g-&gt;<a class="code" href="structtw__rng__stream.html#ae50417100717d3cd46dacb1d7a01e6bf">Cg</a>[1] = s;</div>
<div class="line"><a name="l00415"></a><span class="lineno">  415</span>&#160;  u = u - 4.65661310075985993e-10 * s;</div>
<div class="line"><a name="l00416"></a><span class="lineno">  416</span>&#160;  <span class="keywordflow">if</span>(u &lt; 0)</div>
<div class="line"><a name="l00417"></a><span class="lineno">  417</span>&#160;    u = u + 1.0;</div>
<div class="line"><a name="l00418"></a><span class="lineno">  418</span>&#160;</div>
<div class="line"><a name="l00419"></a><span class="lineno">  419</span>&#160;  s = g-&gt;<a class="code" href="structtw__rng__stream.html#ae50417100717d3cd46dacb1d7a01e6bf">Cg</a>[2];</div>
<div class="line"><a name="l00420"></a><span class="lineno">  420</span>&#160;  k = s / 15499;</div>
<div class="line"><a name="l00421"></a><span class="lineno">  421</span>&#160;  s = 138556 *(s - k * 15499) - k * 3979;</div>
<div class="line"><a name="l00422"></a><span class="lineno">  422</span>&#160;  <span class="keywordflow">if</span>(s &lt; 0.0)</div>
<div class="line"><a name="l00423"></a><span class="lineno">  423</span>&#160;    s = s + 2147483423;</div>
<div class="line"><a name="l00424"></a><span class="lineno">  424</span>&#160;  g-&gt;<a class="code" href="structtw__rng__stream.html#ae50417100717d3cd46dacb1d7a01e6bf">Cg</a>[2] = s;</div>
<div class="line"><a name="l00425"></a><span class="lineno">  425</span>&#160;  u = u + 4.65661336096842131e-10 * s;</div>
<div class="line"><a name="l00426"></a><span class="lineno">  426</span>&#160;  <span class="keywordflow">if</span>(u &gt;= 1.0)</div>
<div class="line"><a name="l00427"></a><span class="lineno">  427</span>&#160;    u = u - 1.0;</div>
<div class="line"><a name="l00428"></a><span class="lineno">  428</span>&#160;</div>
<div class="line"><a name="l00429"></a><span class="lineno">  429</span>&#160;  s = g-&gt;<a class="code" href="structtw__rng__stream.html#ae50417100717d3cd46dacb1d7a01e6bf">Cg</a>[3];</div>
<div class="line"><a name="l00430"></a><span class="lineno">  430</span>&#160;  k = s / 43218;</div>
<div class="line"><a name="l00431"></a><span class="lineno">  431</span>&#160;  s = 49689 *(s - k * 43218) - k * 24121;</div>
<div class="line"><a name="l00432"></a><span class="lineno">  432</span>&#160;  <span class="keywordflow">if</span>(s &lt; 0)</div>
<div class="line"><a name="l00433"></a><span class="lineno">  433</span>&#160;    s = s + 2147483323;</div>
<div class="line"><a name="l00434"></a><span class="lineno">  434</span>&#160;  g-&gt;<a class="code" href="structtw__rng__stream.html#ae50417100717d3cd46dacb1d7a01e6bf">Cg</a>[3] = s;</div>
<div class="line"><a name="l00435"></a><span class="lineno">  435</span>&#160;  u = u - 4.65661357780891134e-10 * s;</div>
<div class="line"><a name="l00436"></a><span class="lineno">  436</span>&#160;  <span class="keywordflow">if</span>(u &lt; 0)</div>
<div class="line"><a name="l00437"></a><span class="lineno">  437</span>&#160;    u = u + 1.0;</div>
<div class="line"><a name="l00438"></a><span class="lineno">  438</span>&#160;</div>
<div class="line"><a name="l00439"></a><span class="lineno">  439</span>&#160;  <span class="keywordflow">return</span> u;</div>
<div class="line"><a name="l00440"></a><span class="lineno">  440</span>&#160;}</div>
<div class="line"><a name="l00441"></a><span class="lineno">  441</span>&#160;</div>
<div class="line"><a name="l00442"></a><span class="lineno">  442</span>&#160;<span class="comment">/*</span></div>
<div class="line"><a name="l00443"></a><span class="lineno">  443</span>&#160;<span class="comment"> * rng_gen_reverse_val</span></div>
<div class="line"><a name="l00444"></a><span class="lineno">  444</span>&#160;<span class="comment"> *</span></div>
<div class="line"><a name="l00445"></a><span class="lineno">  445</span>&#160;<span class="comment"> * computes the reverse sequence, however does not </span></div>
<div class="line"><a name="l00446"></a><span class="lineno">  446</span>&#160;<span class="comment"> * return the uniform value computed as a performance </span></div>
<div class="line"><a name="l00447"></a><span class="lineno">  447</span>&#160;<span class="comment"> * optimization -- Chris Carothers 5/15/98            </span></div>
<div class="line"><a name="l00448"></a><span class="lineno">  448</span>&#160;<span class="comment"> */</span></div>
<div class="line"><a name="l00449"></a><span class="lineno">  449</span>&#160;</div>
<div class="line"><a name="l00450"></a><span class="lineno">  450</span>&#160;<span class="keywordtype">double</span></div>
<div class="line"><a name="l00451"></a><span class="lineno"><a class="line" href="rand-clcg4_8h.html#aa596cdb21f9c04984cc4806d7357d542">  451</a></span>&#160;<a class="code" href="rand-clcg4_8c.html#aa596cdb21f9c04984cc4806d7357d542">rng_gen_reverse_val</a>(<a class="code" href="structtw__rng__stream.html">tw_rng_stream</a> * g)</div>
<div class="line"><a name="l00452"></a><span class="lineno">  452</span>&#160;{</div>
<div class="line"><a name="l00453"></a><span class="lineno">  453</span>&#160;  <span class="keywordtype">long</span> <span class="keywordtype">long</span> *b = rng-&gt;<a class="code" href="structtw__rng.html#a3966a408c9f22a22d2e7ebc387156086">b</a>;</div>
<div class="line"><a name="l00454"></a><span class="lineno">  454</span>&#160;  int32_t *m = rng-&gt;<a class="code" href="structtw__rng.html#ab6cd71dbaaf9c25885b5647c4e1f9d02">m</a>;</div>
<div class="line"><a name="l00455"></a><span class="lineno">  455</span>&#160;  int32_t s;</div>
<div class="line"><a name="l00456"></a><span class="lineno">  456</span>&#160;  <span class="keywordtype">double</span> u;</div>
<div class="line"><a name="l00457"></a><span class="lineno">  457</span>&#160;</div>
<div class="line"><a name="l00458"></a><span class="lineno">  458</span>&#160;  g-&gt;<a class="code" href="structtw__rng__stream.html#ae373cfc091d1ac7361776979d914ad59">count</a>--;</div>
<div class="line"><a name="l00459"></a><span class="lineno">  459</span>&#160;</div>
<div class="line"><a name="l00460"></a><span class="lineno">  460</span>&#160;  u = 0.0;</div>
<div class="line"><a name="l00461"></a><span class="lineno">  461</span>&#160;</div>
<div class="line"><a name="l00462"></a><span class="lineno">  462</span>&#160;  <span class="keywordflow">if</span>(b[0] == 0)</div>
<div class="line"><a name="l00463"></a><span class="lineno">  463</span>&#160;    <a class="code" href="ross-extern_8h.html#a49ed2388aaae26e43280b7909c834aaa">tw_error</a>(<a class="code" href="ross-extern_8h.html#afffdf140f2c492f61546ad783833f127">TW_LOC</a>, <span class="stringliteral">&quot;b values not calculated \n&quot;</span>);</div>
<div class="line"><a name="l00464"></a><span class="lineno">  464</span>&#160;</div>
<div class="line"><a name="l00465"></a><span class="lineno">  465</span>&#160;  s = g-&gt;<a class="code" href="structtw__rng__stream.html#ae50417100717d3cd46dacb1d7a01e6bf">Cg</a>[0];</div>
<div class="line"><a name="l00466"></a><span class="lineno">  466</span>&#160;  s =(b[0] * s) % m[0];</div>
<div class="line"><a name="l00467"></a><span class="lineno">  467</span>&#160;  g-&gt;<a class="code" href="structtw__rng__stream.html#ae50417100717d3cd46dacb1d7a01e6bf">Cg</a>[0] = s;</div>
<div class="line"><a name="l00468"></a><span class="lineno">  468</span>&#160;  u = u + 4.65661287524579692e-10 * s;</div>
<div class="line"><a name="l00469"></a><span class="lineno">  469</span>&#160;</div>
<div class="line"><a name="l00470"></a><span class="lineno">  470</span>&#160;  s = g-&gt;<a class="code" href="structtw__rng__stream.html#ae50417100717d3cd46dacb1d7a01e6bf">Cg</a>[1];</div>
<div class="line"><a name="l00471"></a><span class="lineno">  471</span>&#160;  s =(b[1] * s) % m[1];</div>
<div class="line"><a name="l00472"></a><span class="lineno">  472</span>&#160;  g-&gt;<a class="code" href="structtw__rng__stream.html#ae50417100717d3cd46dacb1d7a01e6bf">Cg</a>[1] = s;</div>
<div class="line"><a name="l00473"></a><span class="lineno">  473</span>&#160;  u = u - 4.65661310075985993e-10 * s;</div>
<div class="line"><a name="l00474"></a><span class="lineno">  474</span>&#160;  <span class="keywordflow">if</span>(u &lt; 0)</div>
<div class="line"><a name="l00475"></a><span class="lineno">  475</span>&#160;    u = u + 1.0;</div>
<div class="line"><a name="l00476"></a><span class="lineno">  476</span>&#160;</div>
<div class="line"><a name="l00477"></a><span class="lineno">  477</span>&#160;  s = g-&gt;<a class="code" href="structtw__rng__stream.html#ae50417100717d3cd46dacb1d7a01e6bf">Cg</a>[2];</div>
<div class="line"><a name="l00478"></a><span class="lineno">  478</span>&#160;  s =(b[2] * s) % m[2];</div>
<div class="line"><a name="l00479"></a><span class="lineno">  479</span>&#160;  g-&gt;<a class="code" href="structtw__rng__stream.html#ae50417100717d3cd46dacb1d7a01e6bf">Cg</a>[2] = s;</div>
<div class="line"><a name="l00480"></a><span class="lineno">  480</span>&#160;  u = u + 4.65661336096842131e-10 * s;</div>
<div class="line"><a name="l00481"></a><span class="lineno">  481</span>&#160;  <span class="keywordflow">if</span>(u &gt;= 1.0)</div>
<div class="line"><a name="l00482"></a><span class="lineno">  482</span>&#160;    u = u - 1.0;</div>
<div class="line"><a name="l00483"></a><span class="lineno">  483</span>&#160;</div>
<div class="line"><a name="l00484"></a><span class="lineno">  484</span>&#160;  s = g-&gt;<a class="code" href="structtw__rng__stream.html#ae50417100717d3cd46dacb1d7a01e6bf">Cg</a>[3];</div>
<div class="line"><a name="l00485"></a><span class="lineno">  485</span>&#160;  s =(b[3] * s) % m[3];</div>
<div class="line"><a name="l00486"></a><span class="lineno">  486</span>&#160;  g-&gt;<a class="code" href="structtw__rng__stream.html#ae50417100717d3cd46dacb1d7a01e6bf">Cg</a>[3] = s;</div>
<div class="line"><a name="l00487"></a><span class="lineno">  487</span>&#160;  u = u - 4.65661357780891134e-10 * s;</div>
<div class="line"><a name="l00488"></a><span class="lineno">  488</span>&#160;  <span class="keywordflow">if</span>(u &lt; 0)</div>
<div class="line"><a name="l00489"></a><span class="lineno">  489</span>&#160;    u = u + 1.0;</div>
<div class="line"><a name="l00490"></a><span class="lineno">  490</span>&#160;</div>
<div class="line"><a name="l00491"></a><span class="lineno">  491</span>&#160;  <span class="keywordflow">return</span> u;</div>
<div class="line"><a name="l00492"></a><span class="lineno">  492</span>&#160;}</div>
<div class="ttc" id="rand-clcg4_8c_html_af6fe13cfedc4e727feee7b28869e036e"><div class="ttname"><a href="rand-clcg4_8c.html#af6fe13cfedc4e727feee7b28869e036e">clamp_seed</a></div><div class="ttdeci">void clamp_seed(uint32_t s[4])</div><div class="ttdef"><b>Definition:</b> <a href="rand-clcg4_8c_source.html#l00153">rand-clcg4.c:153</a></div></div>
<div class="ttc" id="ross-extern_8h_html_afffdf140f2c492f61546ad783833f127"><div class="ttname"><a href="ross-extern_8h.html#afffdf140f2c492f61546ad783833f127">TW_LOC</a></div><div class="ttdeci">#define TW_LOC</div><div class="ttdef"><b>Definition:</b> <a href="ross-extern_8h_source.html#l00164">ross-extern.h:164</a></div></div>
<div class="ttc" id="rand-clcg4_8h_html_a3481da76773e3f2621af85e5865484a5a6f0f1982e69fd30c91d7d3b7097f8912"><div class="ttname"><a href="rand-clcg4_8h.html#a3481da76773e3f2621af85e5865484a5a6f0f1982e69fd30c91d7d3b7097f8912">NewSeed</a></div><div class="ttdef"><b>Definition:</b> <a href="rand-clcg4_8h_source.html#l00027">rand-clcg4.h:27</a></div></div>
<div class="ttc" id="rand-clcg4_8c_html_a3c5ea1f9125b9534b62a80dbfa20dafb"><div class="ttname"><a href="rand-clcg4_8c.html#a3c5ea1f9125b9534b62a80dbfa20dafb">rng_init</a></div><div class="ttdeci">tw_rng * rng_init(int v, int w)</div><div class="ttdef"><b>Definition:</b> <a href="rand-clcg4_8c_source.html#l00338">rand-clcg4.c:338</a></div></div>
<div class="ttc" id="structtw__rng_html"><div class="ttname"><a href="structtw__rng.html">tw_rng</a></div><div class="ttdef"><b>Definition:</b> <a href="rand-clcg4_8h_source.html#l00006">rand-clcg4.h:6</a></div></div>
<div class="ttc" id="ross-extern_8h_html_a49ed2388aaae26e43280b7909c834aaa"><div class="ttname"><a href="ross-extern_8h.html#a49ed2388aaae26e43280b7909c834aaa">tw_error</a></div><div class="ttdeci">void tw_error(const char *file, int line, const char *fmt,...) NORETURN</div><div class="ttdef"><b>Definition:</b> <a href="tw-util_8c_source.html#l00074">tw-util.c:74</a></div></div>
<div class="ttc" id="rand-clcg4_8h_html_a3481da76773e3f2621af85e5865484a5"><div class="ttname"><a href="rand-clcg4_8h.html#a3481da76773e3f2621af85e5865484a5">SeedType</a></div><div class="ttdeci">SeedType</div><div class="ttdef"><b>Definition:</b> <a href="rand-clcg4_8h_source.html#l00025">rand-clcg4.h:25</a></div></div>
<div class="ttc" id="structtw__rng_html_aff3f20edc0b90fd3191487c21b2de697"><div class="ttname"><a href="structtw__rng.html#aff3f20edc0b90fd3191487c21b2de697">tw_rng::a</a></div><div class="ttdeci">int32_t a[4]</div><div class="ttdef"><b>Definition:</b> <a href="rand-clcg4_8h_source.html#l00017">rand-clcg4.h:17</a></div></div>
<div class="ttc" id="structtw__rng_html_a4f5b2d7e79fcebe5c9731b1cf3ee98db"><div class="ttname"><a href="structtw__rng.html#a4f5b2d7e79fcebe5c9731b1cf3ee98db">tw_rng::avw</a></div><div class="ttdeci">int32_t avw[4]</div><div class="ttdef"><b>Definition:</b> <a href="rand-clcg4_8h_source.html#l00019">rand-clcg4.h:19</a></div></div>
<div class="ttc" id="rand-clcg4_8c_html_ab89ab01c93d79746cff9a7ed763b24d3"><div class="ttname"><a href="rand-clcg4_8c.html#ab89ab01c93d79746cff9a7ed763b24d3">MultModM</a></div><div class="ttdeci">int32_t MultModM(int32_t s, int32_t t, int32_t M)</div><div class="ttdoc">Returns(s*t) MOD M. </div><div class="ttdef"><b>Definition:</b> <a href="rand-clcg4_8c_source.html#l00077">rand-clcg4.c:77</a></div></div>
<div class="ttc" id="rand-clcg4_8c_html_adc7b882f7220d319c89b7500c5d2395f"><div class="ttname"><a href="rand-clcg4_8c.html#adc7b882f7220d319c89b7500c5d2395f">rng_gen_val</a></div><div class="ttdeci">double rng_gen_val(tw_rng_stream *g)</div><div class="ttdef"><b>Definition:</b> <a href="rand-clcg4_8c_source.html#l00392">rand-clcg4.c:392</a></div></div>
<div class="ttc" id="rand-clcg4_8c_html_ad0fa533d1469c99034b45cd302a5134c"><div class="ttname"><a href="rand-clcg4_8c.html#ad0fa533d1469c99034b45cd302a5134c">rng_write_state</a></div><div class="ttdeci">void rng_write_state(tw_rng_stream *g, FILE *f)</div><div class="ttdef"><b>Definition:</b> <a href="rand-clcg4_8c_source.html#l00198">rand-clcg4.c:198</a></div></div>
<div class="ttc" id="rand-clcg4_8c_html_a9179f0148c620729309111ad61a3fa5b"><div class="ttname"><a href="rand-clcg4_8c.html#a9179f0148c620729309111ad61a3fa5b">tw_rand_initial_seed</a></div><div class="ttdeci">void tw_rand_initial_seed(tw_rng_stream *g, tw_lpid id)</div><div class="ttdef"><b>Definition:</b> <a href="rand-clcg4_8c_source.html#l00265">rand-clcg4.c:265</a></div></div>
<div class="ttc" id="rand-clcg4_8c_html_ac708ebfb62ad4528bedaeffbe28a3e06"><div class="ttname"><a href="rand-clcg4_8c.html#ac708ebfb62ad4528bedaeffbe28a3e06">S3_MAX</a></div><div class="ttdeci">#define S3_MAX</div><div class="ttdef"><b>Definition:</b> <a href="rand-clcg4_8c_source.html#l00147">rand-clcg4.c:147</a></div></div>
<div class="ttc" id="structtw__rng_html_ad64150f8992ffc7cff059682f5418378"><div class="ttname"><a href="structtw__rng.html#ad64150f8992ffc7cff059682f5418378">tw_rng::aw</a></div><div class="ttdeci">int32_t aw[4]</div><div class="ttdef"><b>Definition:</b> <a href="rand-clcg4_8h_source.html#l00018">rand-clcg4.h:18</a></div></div>
<div class="ttc" id="rand-clcg4_8c_html_a5bb18b1128fadea0f468ea648f656ce1"><div class="ttname"><a href="rand-clcg4_8c.html#a5bb18b1128fadea0f468ea648f656ce1">rng</a></div><div class="ttdeci">static tw_rng * rng</div><div class="ttdef"><b>Definition:</b> <a href="rand-clcg4_8c_source.html#l00023">rand-clcg4.c:23</a></div></div>
<div class="ttc" id="rand-clcg4_8c_html_a2888609ddc8a467786e86846bec7ddb9"><div class="ttname"><a href="rand-clcg4_8c.html#a2888609ddc8a467786e86846bec7ddb9">S0_MAX</a></div><div class="ttdeci">#define S0_MAX</div><div class="ttdef"><b>Definition:</b> <a href="rand-clcg4_8c_source.html#l00144">rand-clcg4.c:144</a></div></div>
<div class="ttc" id="ross_8h_html_a911515dff2fb79886a16ed44df5f9d20"><div class="ttname"><a href="ross_8h.html#a911515dff2fb79886a16ed44df5f9d20">tw_lpid</a></div><div class="ttdeci">uint64_t tw_lpid</div><div class="ttdef"><b>Definition:</b> <a href="ross_8h_source.html#l00160">ross.h:160</a></div></div>
<div class="ttc" id="rand-clcg4_8c_html_a0b3e49b839fa2dd23a2fddf31aaf6ad2"><div class="ttname"><a href="rand-clcg4_8c.html#a0b3e49b839fa2dd23a2fddf31aaf6ad2">S1_MAX</a></div><div class="ttdeci">#define S1_MAX</div><div class="ttdef"><b>Definition:</b> <a href="rand-clcg4_8c_source.html#l00145">rand-clcg4.c:145</a></div></div>
<div class="ttc" id="rand-clcg4_8c_html_abec92cc72a096640b821b8cd56a02495"><div class="ttname"><a href="rand-clcg4_8c.html#abec92cc72a096640b821b8cd56a02495">H</a></div><div class="ttdeci">#define H</div><div class="ttdef"><b>Definition:</b> <a href="rand-clcg4_8c_source.html#l00020">rand-clcg4.c:20</a></div></div>
<div class="ttc" id="structtw__lp_html_a3e75d9bf7796c16870ad686564ee4c20"><div class="ttname"><a href="structtw__lp.html#a3e75d9bf7796c16870ad686564ee4c20">tw_lp::gid</a></div><div class="ttdeci">tw_lpid gid</div><div class="ttdoc">global LP id </div><div class="ttdef"><b>Definition:</b> <a href="ross-types_8h_source.html#l00306">ross-types.h:306</a></div></div>
<div class="ttc" id="ross-extern_8h_html_a11d3d6ea7b82645df5fc7011a407c355"><div class="ttname"><a href="ross-extern_8h.html#a11d3d6ea7b82645df5fc7011a407c355">g_tw_rng_seed</a></div><div class="ttdeci">tw_seed g_tw_rng_seed</div><div class="ttdef"><b>Definition:</b> <a href="ross-global_8c_source.html#l00031">ross-global.c:31</a></div></div>
<div class="ttc" id="structtw__rng_html_a3966a408c9f22a22d2e7ebc387156086"><div class="ttname"><a href="structtw__rng.html#a3966a408c9f22a22d2e7ebc387156086">tw_rng::b</a></div><div class="ttdeci">long long b[4]</div><div class="ttdef"><b>Definition:</b> <a href="rand-clcg4_8h_source.html#l00011">rand-clcg4.h:11</a></div></div>
<div class="ttc" id="rand-clcg4_8c_html_a5a44a44cb9212f0097f9beec79cc8d41"><div class="ttname"><a href="rand-clcg4_8c.html#a5a44a44cb9212f0097f9beec79cc8d41">S2_MAX</a></div><div class="ttdeci">#define S2_MAX</div><div class="ttdef"><b>Definition:</b> <a href="rand-clcg4_8c_source.html#l00146">rand-clcg4.c:146</a></div></div>
<div class="ttc" id="rand-clcg4_8c_html_aaf057a5d3e55f2174f5100a2058f37c5"><div class="ttname"><a href="rand-clcg4_8c.html#aaf057a5d3e55f2174f5100a2058f37c5">rng_get_state</a></div><div class="ttdeci">void rng_get_state(tw_rng_stream *g, uint32_t s[4])</div><div class="ttdef"><b>Definition:</b> <a href="rand-clcg4_8c_source.html#l00212">rand-clcg4.c:212</a></div></div>
<div class="ttc" id="structtw__rng_html_ada588d3e377c9f0d10fb22406f5b47c6"><div class="ttname"><a href="structtw__rng.html#ada588d3e377c9f0d10fb22406f5b47c6">tw_rng::seed</a></div><div class="ttdeci">int32_t seed[4]</div><div class="ttdef"><b>Definition:</b> <a href="rand-clcg4_8h_source.html#l00022">rand-clcg4.h:22</a></div></div>
<div class="ttc" id="structtw__rng__stream_html_aa90b1f32cf686365db9070d700ca0a38"><div class="ttname"><a href="structtw__rng__stream.html#aa90b1f32cf686365db9070d700ca0a38">tw_rng_stream::Ig</a></div><div class="ttdeci">int32_t Ig[4]</div><div class="ttdef"><b>Definition:</b> <a href="rand-clcg4_8h_source.html#l00035">rand-clcg4.h:35</a></div></div>
<div class="ttc" id="ross_8h_html"><div class="ttname"><a href="ross_8h.html">ross.h</a></div></div>
<div class="ttc" id="rand-clcg4_8h_html_a3481da76773e3f2621af85e5865484a5a36cbdbd62f02fd7688d6dfc56c5b9670"><div class="ttname"><a href="rand-clcg4_8h.html#a3481da76773e3f2621af85e5865484a5a36cbdbd62f02fd7688d6dfc56c5b9670">InitialSeed</a></div><div class="ttdef"><b>Definition:</b> <a href="rand-clcg4_8h_source.html#l00027">rand-clcg4.h:27</a></div></div>
<div class="ttc" id="rand-clcg4_8h_html_a3481da76773e3f2621af85e5865484a5a87c13b1e6a443e099ca4aa828db9e7a3"><div class="ttname"><a href="rand-clcg4_8h.html#a3481da76773e3f2621af85e5865484a5a87c13b1e6a443e099ca4aa828db9e7a3">LastSeed</a></div><div class="ttdef"><b>Definition:</b> <a href="rand-clcg4_8h_source.html#l00027">rand-clcg4.h:27</a></div></div>
<div class="ttc" id="ross-extern_8h_html_a99999f4e30648f5007c61723c95a2d0e"><div class="ttname"><a href="ross-extern_8h.html#a99999f4e30648f5007c61723c95a2d0e">g_tw_nRNG_per_lp</a></div><div class="ttdeci">unsigned int g_tw_nRNG_per_lp</div><div class="ttdef"><b>Definition:</b> <a href="ross-global_8c_source.html#l00029">ross-global.c:29</a></div></div>
<div class="ttc" id="structtw__rng__stream_html_a61563015a585f69852a131cb62a4a5b1"><div class="ttname"><a href="structtw__rng__stream.html#a61563015a585f69852a131cb62a4a5b1">tw_rng_stream::Lg</a></div><div class="ttdeci">int32_t Lg[4]</div><div class="ttdef"><b>Definition:</b> <a href="rand-clcg4_8h_source.html#l00036">rand-clcg4.h:36</a></div></div>
<div class="ttc" id="rand-clcg4_8c_html_a29a7a9b6a079ced3ba8339a66856fa82"><div class="ttname"><a href="rand-clcg4_8c.html#a29a7a9b6a079ced3ba8339a66856fa82">rng_init_generator</a></div><div class="ttdeci">void rng_init_generator(tw_rng_stream *g, SeedType Where)</div><div class="ttdef"><b>Definition:</b> <a href="rand-clcg4_8c_source.html#l00239">rand-clcg4.c:239</a></div></div>
<div class="ttc" id="rand-clcg4_8c_html_a9c251d8e1876abc122b55fe4156ce14c"><div class="ttname"><a href="rand-clcg4_8c.html#a9c251d8e1876abc122b55fe4156ce14c">rng_set_seed</a></div><div class="ttdeci">void rng_set_seed(tw_rng_stream *g, uint32_t s[4])</div><div class="ttdef"><b>Definition:</b> <a href="rand-clcg4_8c_source.html#l00181">rand-clcg4.c:181</a></div></div>
<div class="ttc" id="rand-clcg4_8c_html_aa97731a64b1cf687109f7a5f8f519c68"><div class="ttname"><a href="rand-clcg4_8c.html#aa97731a64b1cf687109f7a5f8f519c68">rng_put_state</a></div><div class="ttdeci">void rng_put_state(tw_rng_stream *g, uint32_t s[4])</div><div class="ttdef"><b>Definition:</b> <a href="rand-clcg4_8c_source.html#l00226">rand-clcg4.c:226</a></div></div>
<div class="ttc" id="rand-clcg4_8c_html_aa596cdb21f9c04984cc4806d7357d542"><div class="ttname"><a href="rand-clcg4_8c.html#aa596cdb21f9c04984cc4806d7357d542">rng_gen_reverse_val</a></div><div class="ttdeci">double rng_gen_reverse_val(tw_rng_stream *g)</div><div class="ttdef"><b>Definition:</b> <a href="rand-clcg4_8c_source.html#l00451">rand-clcg4.c:451</a></div></div>
<div class="ttc" id="structtw__rng__stream_html_ae373cfc091d1ac7361776979d914ad59"><div class="ttname"><a href="structtw__rng__stream.html#ae373cfc091d1ac7361776979d914ad59">tw_rng_stream::count</a></div><div class="ttdeci">unsigned long count</div><div class="ttdef"><b>Definition:</b> <a href="rand-clcg4_8h_source.html#l00034">rand-clcg4.h:34</a></div></div>
<div class="ttc" id="structtw__lp_html_af7e48c8a55f5fff35746ac044df1902f"><div class="ttname"><a href="structtw__lp.html#af7e48c8a55f5fff35746ac044df1902f">tw_lp::rng</a></div><div class="ttdeci">tw_rng_stream * rng</div><div class="ttdoc">RNG stream array for this LP. </div><div class="ttdef"><b>Definition:</b> <a href="ross-types_8h_source.html#l00317">ross-types.h:317</a></div></div>
<div class="ttc" id="structtw__rng_html_ab6cd71dbaaf9c25885b5647c4e1f9d02"><div class="ttname"><a href="structtw__rng.html#ab6cd71dbaaf9c25885b5647c4e1f9d02">tw_rng::m</a></div><div class="ttdeci">int32_t m[4]</div><div class="ttdef"><b>Definition:</b> <a href="rand-clcg4_8h_source.html#l00016">rand-clcg4.h:16</a></div></div>
<div class="ttc" id="structtw__rng__stream_html_ae50417100717d3cd46dacb1d7a01e6bf"><div class="ttname"><a href="structtw__rng__stream.html#ae50417100717d3cd46dacb1d7a01e6bf">tw_rng_stream::Cg</a></div><div class="ttdeci">int32_t Cg[4]</div><div class="ttdef"><b>Definition:</b> <a href="rand-clcg4_8h_source.html#l00037">rand-clcg4.h:37</a></div></div>
<div class="ttc" id="rand-clcg4_8c_html_a96c2fb54945c414a31ce0512c39622b0"><div class="ttname"><a href="rand-clcg4_8c.html#a96c2fb54945c414a31ce0512c39622b0">tw_rand_init_streams</a></div><div class="ttdeci">void tw_rand_init_streams(tw_lp *lp, unsigned int nstreams)</div><div class="ttdef"><b>Definition:</b> <a href="rand-clcg4_8c_source.html#l00320">rand-clcg4.c:320</a></div></div>
<div class="ttc" id="rand-clcg4_8c_html_ac243c54261ecf17a1b848e9b9a2025a9"><div class="ttname"><a href="rand-clcg4_8c.html#ac243c54261ecf17a1b848e9b9a2025a9">FindB</a></div><div class="ttdeci">long long FindB(long long a, long long k, long long m)</div><div class="ttdoc">Find B to run CLCG4 backwards. </div><div class="ttdef"><b>Definition:</b> <a href="rand-clcg4_8c_source.html#l00037">rand-clcg4.c:37</a></div></div>
<div class="ttc" id="rand-clcg4_8c_html_ac6110bfd61d392bab210049c5f6f8fa3"><div class="ttname"><a href="rand-clcg4_8c.html#ac6110bfd61d392bab210049c5f6f8fa3">seed</a></div><div class="ttdeci">int32_t seed[4]</div><div class="ttdef"><b>Definition:</b> <a href="rand-clcg4_8c_source.html#l00026">rand-clcg4.c:26</a></div></div>
<div class="ttc" id="ross-extern_8h_html_a445b1ed30cd2e8ef1f837e782a339a70"><div class="ttname"><a href="ross-extern_8h.html#a445b1ed30cd2e8ef1f837e782a339a70">tw_calloc</a></div><div class="ttdeci">void * tw_calloc(const char *file, int line, const char *for_who, size_t e_sz, size_t n)</div><div class="ttdef"><b>Definition:</b> <a href="tw-util_8c_source.html#l00203">tw-util.c:203</a></div></div>
<div class="ttc" id="structtw__rng__stream_html"><div class="ttname"><a href="structtw__rng__stream.html">tw_rng_stream</a></div><div class="ttdef"><b>Definition:</b> <a href="rand-clcg4_8h_source.html#l00032">rand-clcg4.h:32</a></div></div>
<div class="ttc" id="structtw__lp_html"><div class="ttname"><a href="structtw__lp.html">tw_lp</a></div><div class="ttdoc">LP State Structure. </div><div class="ttdef"><b>Definition:</b> <a href="ross-types_8h_source.html#l00304">ross-types.h:304</a></div></div>
</div><!-- fragment --></div><!-- contents -->
<!-- start footer part -->
<hr class="footer"/><address class="footer"><small>
Generated by &#160;<a href="http://www.doxygen.org/index.html">
<img class="footer" src="doxygen.png" alt="doxygen"/>
</a> 1.8.6
</small></address>
</body>
</html>
